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Abstract. We investigate the origin of the abundance 
variations observed among similarly aged F and G dwarfs 
in the local Galactic disk. We argue that orbital diffusion 
of stars in combination with radial abundance gradients 
is probably insufficient to explain these variations. 

We show that episodic and local infall of metal- 
deficient gas can provide an adequate explanation for iron 
and oxygen variations as large as A[M/H] ~0.6 dex among 
stars formed at a given age in the solar neighbourhood 
(SNBH). However, such models appear inconsistent with 
the observations because they: 1) result in current disk 
ISM abundances that are too high compared to the obser- 
vations, 2) predict stellar abundance variations to increase 
with the lifetime of the disk, and 3) do not show substan- 
tial scatter in the [Fe/H] vs. [O/H] relation. Notwithstand- 
ing, our results do suggest that metal-deficient gas infall 
plays an important role in regulating the chemical evolu- 
tion of the Galactic disk. 

We demonstrate that sequential enrichment by succes- 
sive stellar generations within individual gas clouds can 
account for substantial abundance variations as well. How- 
ever, such models are inconsistent with the observations 
because they: 1) are unable to account for the full magni- 
tude of the observed variations, in particular for [Fe/H], 2) 
predict stellar abundance variations to decrease with the 
lifetime of the disk, and 3) result in current abundances 
far below the typical abundances observed in the local disk 
ISM. 

We present arguments in support of combined infall of 
metal-deficient gas and sequential enrichment by succes- 
sive stellar generations in the local Galactic disk ISM. We 
show that galactic chemical evolution models which take 
into account these processes simultaneously are consis- 
tent with both the observed abundance variations among 
similarly aged F and G dwarfs in the SNBH and the 
abundances observed in the local disk ISM. For reason- 
able choices of parameters, these models can reproduce 



A[M/H] for individual elements M = C, O, Fe, Mg, Al, and 
Si as well as the scatter observed in abundance-abundance 
relations like [O/Fe]. For the same models, the contri- 
bution of sequential stellar enrichment to the magnitude 
of the observed abundance variations can be as large as 
-50%. 

We discuss the impact of sequential stellar enrichment 
and episodic infall of metal-deficient gas on the inhomo- 
geneous chemical evolution of the Galactic disk. 

Key words: Galaxy: chemical evolution, abundances, so- 
lar neighbourhood - ISM: abundances - Galaxies: ISM 
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1. Introduction 

The chemical enrichment of the interstellar medium (ISM) 
by successive generations of stars is a key issue in un- 
derstanding the chemical evolution of galaxies in general, 
and the formation history and abundance distributions of 
the stellar populations in our Galaxy in particular. Ob- 
servational studies related to the heavy element enrich- 
ment of the local Galactic disk have long shown that 
stars of similar age exhibit large abundance variations 
(e.g. Mayor 1976; Twarog 1980a; Twarog & Wheeler 1982; 
Carlberg et al. 1985; Gilmore 1989; Klochkova et al. 1989; 
Schuster & Nissen 1989; Meusinger et al. 1991). Recently, 
Edvardsson et al. (1993a) presented accurate abundance 
data for nearly 200 F and G main-sequence dwarfs in the 
solar neighbourhood (SNBH). Their spectroscopic data, 
analysed with up-to-date input physics, confirms abun- 
dance variations as large as ~0.6 dex in A[M/H] (where 
M=Fe,0,Mg,Al,Si) among similarly aged stars. In con- 
trast to previous understanding, these variations are much 
in excess of experimental uncertainties and demonstrate 
that the abundance spread for stars born at roughly the 
same galactocentric distance is similar in magnitude to 
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the overall increase in metallicity during the lifetime of 
the disk. 

Additional support for the existence of large abun- 
dance inhomogeneities in the Galactic disk has been pro- 
vided by studies of stars in open clusters (e.g. Nissen 1988; 
Bocsgaard 1989; Lambert 1989; Garcia-Lopcz et al. 1993; 
Friel & Janes 1993: Carraro & Chiosi 1994) and B stars 
in star forming regions in the SNBH (e.g. Gies & Lambert 
1992; Cunha & Lambert 1992). These studies show that 
the concept of a well-defined tight age-metallicity relation 
(AMR) for the Galactic disk ISM is unfounded (Edmunds 
1993) and that the chemical enrichment of the disk has 
been inhomogeneous on time scales as short as ~10 8 — 10 9 
yr. Similar studies of objects in the Magellanic clouds (e.g. 
Cohen 1982; Da Costa 1991; Olsewski et al. 1991) and 
dwarf galaxies (Pilyugin 1992; Kunth et al. 1994; Thuan 
et al. 1995) suggest that inhomogeneous chemical evolu- 
tion is a common phenomenon in nearby galaxies as well. 

The origin of the abundance variations observed in the 
local Galactic disk is investigated in this paper. Clearly, 
large abundance variations in the ISM on time scales at 
least an order of magnitude shorter than the lifetime of 
the disk cannot be reproduced by simple galactic evolu- 
tion models incorporating monotonously increasing age- 
metallicity relations (AMR). In the past few years, various 
ideas have been put forward as possible explanations for 
the intrinsic abundance variations among similarly aged 
stars: 

• stellar orbital diffusion in combination with radial 
abundance gradients in the Galactic disk (e.g. Fran- 
cois & Matteucci 1993; Wielen, Fuchs, & Dettbarn 
1996); 

• sequential enrichment by successive stellar genera- 
tions (e.g. Edmunds 1975; Olive & Schramm 1982; 
Gilmore 1989; Gilmore & Wyse 1991; Cunha & Lam- 
bert 1992; Roy and Kunth 1995); 

• local infall of metal-poor gas (e.g. Edvardsson et 
al. 1993a; Roy & Kunth 1995; Pilyugin & Edmunds 
1995a); 

• cloud motions in the ISM (Bateman & Larson 1993); 

• inefficient mixing in the disk ISM: isolated chemical 
evolution of individual parcels of interstellar gas dur- 
ing considerable fractions of the lifetime of the disk 
(Lennon et al. 1990; Wilmes & Koppen 1995); 

• major galaxy merger events resulting in multiple stel- 
lar populations in the Galactic disk (Strobel 1991; Pi- 
lyugin & Edmunds 1995b); 

• chemical fractionization processes such as grain for- 
mation (e.g. Henning & Giirtler 1986) and/or ele- 
ment diffusion (e.g. Bahcall & Pinsonneault 1995) so 
that measured abundances do not reflect initial stellar 
abundances. 

As discussed by Edvardsson et al. (1993a) stellar or- 
bital diffusion is probably inadequate as main explanation 
for the observed abundance variations. However, recently 



Wielen et al. (1996) claimed that stars can be born at 
galactocentric distances very different from those derived 
using their present-day orbits. In this case, a major frac- 
tion of the observed abundance scatter could be due to 
stellar orbital diffusion. Since it appears unlikely that dif- 
fusion can explain the observed abundance variations for 
all stars in the Edvardsson et al. sample (see below) as well 
as those observed among young stars (e.g. present in star 
forming regions) other processes are probably important 
as well. 

Based on the assumption of short mixing time scales of 
<~ 10 7 yr in the local disk ISM, we argue that the underly- 
ing physical mechanisms causing the observed abundance 
inhomogeneities and those initiating star formation in the 
disk ISM arc the same. No observational support exists 
for chemical fractionization in low mass F and G main- 
sequence stars. Therefore, the scatter in stellar metallici- 
ties probably reflects the original inhomogeneities in the 
interstellar gas (e.g. Gilmore 1989). We here restrict our- 
selves mainly to the processes of sequential stellar enrich- 
ment and episodic infall of gas onto the Galactic disk as 
possible explanation for the observed stellar abundance 
variations. Since these processes are observed to operate 
simultaneously in the SNBH (see below), it is important 
to investigate their combined effect on the chemical evo- 
lution of the Galactic disk. 

The process of star formation initiated by stars formed 
during a preceding star formation event nearby, is known 
as sequential star formation. Support for sequential star 
formation in the SNBH is provided by observations of spa- 
tially separated subgroups of OB stars that appear aligned 
in a sequence of ages in many OB associations (e.g. Blaauw 
1991) and by obervations of stars forming at the interfaces 
of Hn regions and their surrounding molecular clouds (e.g. 
Genzel & Stutzki 1989; Pismis 1990; Goldsmith 1995). Se- 
quential star formation may be induced by the blast waves 
of nearby supernova explosions compressing the ambient 
ISM (e.g. Ogelman & Maran 1976) and/or by propagat- 
ing ionization and shock fronts from an OB association 
causing the gravitational collapse of a nearby molecular 
cloud (e.g. Elmegreen & Lada 1977). In either case, effi- 
cient self-enrichment through mixing of enriched material 
by successive generations of massive stars is expected. 

On the other hand, stellar abundance variations can 
be attributed to infall of relatively unprocessed gas and 
star formation within the accreted material before effi- 
cient mixing wipes out any local chemical inhomogeneities 
(e.g. Edvardsson et al. 1993a; Pilyugin & Edmunds 1995a). 
Observational support for star formation in the SNBH ini- 
tiated by infall of high velocity clouds, has recently been 
presented by Lepine & Duvert (1994). These authors claim 
that episodic gas infall is a dominant process in the local 
disk ISM and is associated with all prominent star form- 
ing molecular clouds seen near the Sun (see Sect. 5.2). 
Furthermore, ongoing gas infall has been emphasized by 
models of dissipative protogalactic collapse (Larson 1969, 
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1976) and on the basis of time scale arguments of gas con- 
sumption in the local disk (Larson et al. 1980; Kennicutt 
1983). 

In this paper, we present a chemical evolution model 
for a star forming gas cloud which incorporates stellar 
enrichment and mixing processes (including infall) and 
which allows for temporal and/or spatial inhomogeneities 
in the ISM. This study differs from previous work (e.g. 
Pilyugin & Edmunds 1995a) in that we investigate in de- 
tail the combined effect of metal deficient gas infall and 
sequential stellar enrichment by successive stellar genera- 
tions on the chemical evolution of multiple gas clouds in 
the Galactic disk. In particular, each gas cloud is allowed 
to follow its individual star formation, mixing, and in- 
fall history, as is suggested by the observations. With this 
model, we fit the stellar abundance variations and current 
local ISM abundances of C, O, Fe, Mg, Al, and Si ob- 
served in the SNBH, the present gas-to-total mass-ratio, 
and actual star formation and supernova rates. We note 
that previous investigations were restricted to abundance 
variations in oxygen and iron only. 

We will show that models taking into account the 
above processes simultaneously are in good agreement 
with the observations and provide an adequate explana- 
tion for the stellar abundance variations with respect to 
the mean abundances observed in the local disk. Further- 
more, we will argue that the contribution of sequential 
stellar enrichment to the magnitude of the observed stellar 
abundance variations can be as large as ^50%, i.e. much 
larger than suggested by previous investigations (e.g. Pi- 
lyugin & Edmunds 1995a; Wilmes and Koppen 1995). Cor- 
responding theoretical age and abundance-distributions 
related to the G-dwarf problem will be discussed in a sep- 
arate paper. 

The paper is organized as follows. In Sect. 2, we briefly 
review observations related to the inhomogeneous heavy 
clement enrichment of the local Galactic disk ISM. In 
Sect. 3, we describe characteristics of the inhomogeneous 
chemical evolution model proposed for the Galactic disk 
(model equations and details are given in the Appendix 
to the electronic version of this paper). In Sect. 4, we 
present model results for episodic infall of metal-deficient 
gas and sequential stellar enrichment, and examine which 
of these mechanisms can account satisfactorily for the ob- 
servations. In Sect. 5, we discuss these results in the more 
general context of the chemical evolution of the Galactic 
disk and adduce both observational arguments in support 
of sequential star formation and metal deficient gas infall 
in the local disk ISM. 



2. Inhomogeneous chemical evolution of the local 
Galactic disk: observations 



2.1. Main-sequence F and G dwarfs 

We concentrate on the abundance data of nearly 200 main- 
sequence field F and G dwarfs with actual distances ^70 
pc from the Sun as recently presented by Edvardsson et al. 
(1993a; hereafter EDV). This sample provides the largest 
sample of stars available to date for studies related to the 
chemical evolution of the local disk. Fig. 1 displays all F 
and G dwarfs for which both [O/H] and [Fe/H] abundance- 
ratios have been determined by EDV. 

Large abundance variations of ~0.9 dex in [Fe/H] and 
^0.7 dex in [O/H] among stars of a given age are seen to 
be present (abundance variations in e.g. Mg, Al, and Si re- 
semble those in [Fe/H]). At intermediate stellar ages, these 
variations are no doubt significant since typical observa- 
tional errors are ~0.1 dex both in [M/H] and log(Age) 
(see EDV). At ages in excess of ~15 Gyr and less than ^2 
Gyr, the data probably are undersampled (see EDV). Note 
that the sample is biased against old, high-metallicity stars 
through the minimum T e g limit assumed by EDV. 

The observed spread in [Fe/H] is tightly correlated 
with that in [O/H]. This suggests that different nucleo- 
synthesis sites, which contribute different elements to the 
initial abundances in stars, mix their products together 
well. Furthermore, this suggests that stellar abundance 
variations for different elements are due to the same pro- 
cess. Current observations support the idea that the mag- 
nitude of the stellar abundance variations has remained 
constant over the lifetime of the disk (see also Mayor 
1976; Twarog 1980a; Meusinger et al. 1991; Carraro & 
Chiosi 1994). In the following, we will assume that these 
abundance variations are randomly distributed within the 
metallicity range observed at a given stellar age. This is 
particularly important when considering possible expla- 
nations for the observed stellar abundance variations in 
detail (see Sect. 4; cf. Wielen et al. 1996). 



Ages, abundances, and kinematical properties of the 
dwarfs belonging to the EDV sample are consistent with 
earlier investigations (e.g. Twarog 1980a; Carlberg et al. 
1985; Meusinger et al. 1991). An extensive discussion of 
the possible sources of errors in the abundance and age 
analysis as well as several consistency checks can be found 
in Edvardsson et al. (1993a, b). Errors due to data reduc- 
tion uncertainties are estimated to lead to errors of at 
most 0.05 to 0.1 dex in abundance ratios [M/Fe] as well 
as in [Fe/H] (see EDV). These errors are not expected 
to vary in a systematic manner with the derived stellar 
abundances and corresponding corrections will probably 
not reduce the observed variations. Edvardsson et al. esti- 
mated errors in relative ages of ~25% for stars with similar 
abundances (absolute errors may be considerably larger). 
Thus, ages of stars as old as the Sun are estimated to be 
accurate within ~1— 2 Gyr. We conclude that errors in 
the abundances and ages of the sample stars are unlikely 
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Fig. 1. Observed iron, oxygen, and magnesium abundance ratios for main-sequence F and G dwarfs in the solar neighbourhood 
(data from Edvardsson et al. 1993a). Open circles represent stars with mean stellar galactocentric distances at birth within 0.5 
kpc from the Sun (Rq = 8.4 kpc). Full dots indicate stars with average distances within ~2 kpc from the Sun. Typical errors 
are indicated at the bottom right of each panel. Note that the abundances of the most metal-poor disk stars included in this 
sample resemble those of metal-rich halo dwarfs and giants (e.g. Bessell et al. 1991; Gratton & Sneden 1991; Nissen et al. 1994). 
We assumed solar abundance ratios by number of 10 log (O/H) = -3.13, 10 log (Fe/H) Q = -4.51, and 10 log (Mg/H) Q = -4.42 
and a hydrogen mass fraction in the Sun of 0.68 (see Anders & Grevesse 1989; Grevesse & Noels 1993). Top panels: Distributions 
of [Fe/H] (left) and [O/H] abundance ratios vs. galactic age. Bottom panels: [Fe/H] vs. [O/H] (left) and [Mg/H] vs. [O/H] 



to account for the observed abundance variations, at least 
for stars with intermediate ages of ^5 Gyr. 

Knowledge of the formation sites of the sample stars 
is important to decide whether or not orbital diffusion 
in combination with radial abundance gradients in the 
Galactic disk can provide an adequate explanation for the 
observed abundance variations. Galactocentric distances 
of the sample stars at birth were obtained using stellar 
orbits reconstructed from their present-day galactocentric 



distances, proper motions, and radial velocities, and using 
both theoretical and empirical models for the Galactic po- 
tential as discussed by EDV. Accordingly, nearly 85% of 
the sample stars were found to have mean galactocentric 
distances i? m at birth within 1 kpc from the Sun (assum- 
ing Rq ~ 8.4 kpc at present). However, predictions of the 
diffusion of stellar orbits in space, based on the observed 
relation between velocity dispersion and age for nearby 
stars, suggest that many stars may have been formed at 
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galactocentric distances as large as —4 kpc from where 
they are nowadays observed in the SNBH (Wielen et al. 
1996). In cither case, these nearby stars trace the evolu- 
tion of the Galactic disk ISM over a much wider range in 
galactocentric distance then they are observed. 

As an independent test to examine whether stellar or- 
bital diffusion can be the main cause for the observed 
abundance variations, we translated stellar abundance de- 
viations A[M/H] from the mean abundances of similar 
aged stars born at <~ Rq into galactocentric distance dif- 
ferences R m — Rq. This was done independently for [Fe/H] 
and [O/H] abundance ratios assuming present-day local 
radial abundance gradients of — 0.07±0.015 dex kpc^ 1 in 
[O/H] (e.g. Shaver et al. 1983; Grcnon 1987; see also Wil- 
son & Matteucci 1992) and -0.1 dex kpc" 1 in [Fe/H] (see 
e.g. EDV). Clearly, distances R m based on oxygen and iron 
abundances are expected to be similar (e.g. Ai?^' Fc ^ 1 
kpc) when orbital diffusion is important for the observed 
stellar abundance variations. 

In this manner, we find that Ai?° ,Fo <; 1 kpc for 47 
stars in the EDV sample (i.e. —56%). Similarly, we find 
Ai?°< Fo ^ 2 , 3, and 4 kpc, for -31, 14, and 8% of the sam- 
ple stars, respectively. We note that the derived values of 
Ai?° ,Fo are insensitive to the stellar age but depend on 
the assumed radial abundance gradients as well as on the 
mean [M/H] vs. age relations adopted for stars born at 
Rq . Although a detailed investigation of the uncertainties 
involved is beyond the scope of this paper (e.g. the vari- 
ation of abundance gradients with disk age; cf. Grenon 
1987), we estimate that AR^ Fc ^ 0.8 (1.5) kpc for typi- 
cal errors of 0.05 (0.1) dex in both [Fc/H] and [O/H] for 
most of the sample star (assuming a gaussian error dis- 
tribution). This suggests that a substantial part of the 
observed abundance variations is difficult to explain by 
stellar orbital diffusion only. 

Also, Edvardsson et al. argued that the magnitude of 
the observed variations will reduce to A [M/H] — 0.3 dex 
if one accounts properly for systematic errors and possi- 
ble effects of stellar orbital diffusion. However, a reduction 
of the abundance spread among field dwarfs to A [Fe/H] 

— 0.3 dex seems contradicted by the observed variations 
of A [Fe/H] — 0.5 ±0.1 dex among similarly aged open 
clusters after correcting for radial abundance gradients 
across the Galactic plane (Carraro & Chiosi 1994). Apart 
from this, such a reduction appears inconsistent with the 
observed abundance spread of [O/H] ^0.4 dex among B 
stars at a given galactocentric radius between 7 and 16 kpc 
in the disk (Gehren et al. 1985; Kaufer et al. 1994) and 
with the large abundance variations of —0.7 dex observed 
among young open clusters over a distance scale of only 

— 1 kpc at a galactocentric radius of —13 kpc (Rolleston 
et al. 1994). 

What fraction of the current disk stellar population 
actually formed in the Galactic halo depends on the de- 
tailed dynamical evolution of the disk which is not well 
known (e.g. Pagel & Tautvaisiene 1995). However, most 



stars in the EDV sample have derived maximum distances 
from the Galactic plane at birth of h max <0.5 kpc. This 
largely excludes halo stars from the sample and further 
implies that abundance gradients perpendicular to the 
Galactic plane are inadequate as explanation for the ob- 
served abundance variations (e.g. Carney et al. 1990). 

From these arguments, we conclude that orbital diffu- 
sion of stars from elsewhere in the Galactic disk is proba- 
bly insufficient as explanation for the observed variations 
in [Fe/H] and [O/H] among F and G dwarfs in the SNBH. 
This conclusion is consistent with the finding that abun- 
dance variations for subsamples of stars restricted to be 
born within 1 and 0.5 kpc from the Sun, respectively, are 
similar to those for the complete sample (see EDV; cf. Fig. 
1). Therefore, we believe that differential chemical evolu- 
tion and mixing of interstellar gas must be an important 
cause for the large stellar abundance variations observed 
in the SNBH as well. The abundances of the Sun and of 
open clusters in the Galactic disk fit well into this picture, 
as is argued below. 

2.2. Chemical evolution of the solar neighbourhood 

Detailed comparison of abundances within local Hn re- 
gions and the Sun have shown that oxygen (among other 
heavy elements) is underabundant in the Hll regions by 
about 0.15—0.3 dex (e.g. Shaver et al. 1983; Peimbert 
1987; Baldwin et al. 1991; Osterbrock, Tran & Veilleux 
1992). Also, CNO-abundanccs of Hn-regions and B main- 
sequence stars in the Orion nebula were found smaller 
than corresponding abundances in the Sun (Cunha & 
Lambert 1992; Gies & Lambert 1992). The remarkable re- 
sult that the Sun is metal-rich by ~0. 15—0.2 dex in [O/H] 
compared to its surroundings is also supported by observa- 
tions of B stars in nearby associations and young clusters 
(Fitzsimmons et al. 1990), diffuse interstellar clouds (e.g. 
York et al. 1983), and disk planetary nebulae (de Freitas 
Pacheco 1993; Peimbert et al. 1993). Although abundance 
determinations in the SNBH may be biased towards re- 
gions associated with infall of metal-poor gas or suffer 
from heavy clement depletion by dust, the existence of 
many metal-poor regions in the SNBH would be difficult 
to reconcile with efficient mixing in the local disk ISM 
(e.g. Roy & Kunth 1995). 

The above observations are consistent with the Ed- 
vardsson et al. data which suggest that the Sun is metal- 
rich by 0.2-0.25 dex in [O/H] and by 0.25-0.3 dex in 
[Fe/H] compared to the mean abundances of stars which 
formed in the SNBH -4.5 Gyr ago (cf. Fig. 1). These ob- 
servations support the idea that the Sun is metal-rich for 
its age (see also Steigman 1993) and that abundance inho- 
mogeneities in the local disk ISM did exist. The fact that 
the Sun is metal-rich by a factor of —1.5—2 compared to 
nearby regions currently experiencing star formation may 
be explained by self-enrichment of the gas cloud out of 
which the Sun was born (e.g. Gies & Lambert 1992; Pe- 
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imbcrt ct al. 1993). Alternatively, orbital diffusion of the 
Sun may play an important role (Wielen ct al. 1996). We 
will discuss arguments in support of the former possibility 
in Sect. 5.2. 

2.3. Open clusters 

Variations in [Fe/H] among disk open clusters of a given 
age are known to be larger than any possible trend of 
[Fe/H] with age (e.g. Nissen 1988; Boesgaard 1989; Garcia- 
Lopez ct al. 1993; Friel and Janes 1993). Recently, abun- 
dance variations of ^0.5 ±0.1 dex in [Fc/H] among clus- 
ters of a given age after correcting for the radial abundance 
gradient across the Galactic plane have been reported by 
Carraro & Chiosi (1994). The observed abundance varia- 
tions among open clusters appear somewhat smaller (i.e. 
by ~0.2-0.3 dex in [Fe/H]) than those among field F and 
G dwarfs in the EDV sample. However, the magnitude 
of the observed variations suggests that the processes re- 
sponsible for the abundance inhomogencities among field 
stars in the SNBH and among open clusters widespread 
throughout the Galactic disk may well be the same. 

The lack of a tight age-metallicity relationship for open 
clusters in the Galactic disk suggests that the chemical 
enrichment of the disk ISM has been inhomogeneous on 
time scales less than ~10 8 — 10 9 yr, consistent with the 
abundance variations observed for intermediate age F and 
G dwarfs discussed above. 

3. Model characteristics and assumptions 

In the previous section, we have argued that differential 
chemical evolution and mixing of interstellar gas prob- 
ably provides the main explanation for the large abun- 
dance variations observed among similarly aged stars in 
the SNBH. In this case, abundance inhomogeneities in 
the global disk ISM may result from local mixing of 
mctal-dcficicnt material (e.g. infall) and/or local mixing of 
metal-enhanced material (e.g. stellar enrichment). When 
star formation is initiated within the mixed material be- 
fore any abundance fluctuations are wiped out, these in- 
homogeneities can be recorded by long-living stars. 

Efficient mixing by stellar winds and supernova explo- 
sions is generally accepted to occur within ~T0 7 — 10 8 
yr (e.g. Edmunds 1975; Ciotti et al. 1991; Roy & Kunth 
1995). This suggests that the processes responsible for 
the onset of star formation and those causing substantial 
abundance inhomogeneities in the disk ISM are the same. 
We consider this as a strong argument in favour of sequen- 
tial star formation and/or infall induced star formation as 
the main processes responsible for the observed abundance 
variations (ample observational support for the occurence 
of these processes in the local disk ISM are briefly dis- 
cussed in Sect. 5.2). Obviously, the quantitative effect of 
these processes on stellar abundance variations, relative 



to the mean abundances in the local ISM, depends on the 
detailed chemical evolution of the disk ISM. 

3.1. Model description 

We present a model for the inhomogeneous chemical evo- 
lution of a star forming gas cloud. The basis for this model 
forms the individual star formation history and chemical 
evolution of multiple subclouds that mutually exchange 
interstellar material. We here rcstict ourselves to a brief 
outline of the basic assumptions and model characteristics. 
A more detailed description of the equations and input 
physics used is given in the (Appendix to the) electronic 
version of this paper. 

We start from a homogeneous gas cloud with total 
mass M c i. At any time the cloud is subdivided into N sc i 
star forming, active subclouds (with corresponding masses 
M* cl ) and a quiescent, inactive cloud part (with mass M qc i) 
not experiencing star formation. Each subcloud is allowed 
to follow its individual star formation, infall, and mixing 
history. Infall of matter is considered by allowing episodic 
mixing of metal-deficient material to each subcloud sepa- 
rately. 

The adopted set of processes that modify the distri- 
bution of gas and stars within a star forming region of a 
molecular cloud are illustrated in Fig. 2. Different subfig- 
ures refer to the following processes: 

(a) subcloud formation from the inactive cloud ISM 
(and/or from infalling material); 

(b) conversion of gas into stars (star formation at dis- 
tinct subcloud cores); 

(c) ejection of material by stars to their immediate sur- 
roundings; 

(d) mixing of dispersed core material with subcloud af- 
ter star formation event; 

(e) break up of entire subcloud, mixing with inactive 
cloud ISM, and induced star formation; 

(f) enrichment of subcloud by stars not formed within 
the subcloud. 

In our model, the inhomogeneous chemical evolution of a 
star forming gas cloud, consisting of many subclouds, is 
determined by the combined effect of the above processes. 
During a time-interval At, these processes may occur si- 
multaneously within each subcloud. In this manner, the 
initial abundances of a newly formed stellar generation 
are determined by: 1) the enrichment of the subcloud by 
preceding stellar generations, and 2) the mixing history of 
the subcloud with the ambient ISM. 

In brief, the adopted evolution scenario is as follows. 
Subcloud formation (Fig. 2a) is assumed to occur either 
from the inactive cloud ISM and/or from infalling material 
(details related to the infall model will be given in Sect. 
4.3). During the lifetime t cv of the entire system, a total 
number of iV s f star formation events is assumed to occur. 
Each star formation event is assumed to take place in an 



L.B. van den Hoek & T. de Jong: Inhomogeneous chemical evolution of the Galactic disk 



7 






Fig. 2. Schematic model for the inhomogeneous chemical evolution of a star forming cloud: evolutionary sequence of star 
formation, enrichment and mixing processes. Shown is a star forming cloud region. Each of the processes indicated in this 
region may occur in other regions of the cloud as well. Symbols have the following meaning: a subclouds indicated by hatched 
areas, b star formation indicated by asterisks, c stellar enrichment shown as shaded areas enclosing white asterisks, d subcloud 
core dispersal indicated as blanked out area surrounding stars, e break up of entire subcloud and initiation of star formation 
in a nearby subcloud, f arrow indicates stars entering a subcloud from elsewhere. Each of the processes indicated may occur 
frequently during the cloud evolution time t cv 



active subcloud (Fig. 2b) . Each subcloud is allowed to ex- 
perience numerous star formation events and/or to remain 
inactive during a substantial part of its lifetime. Conse- 
quently, each subcloud can be enriched by one or mul- 
tiple star formation events dictating its chemical evolu- 
tion (Fig. 2c). When the active subcloud core is dispersed 
by stellar winds and/or supernova shocks, part of the en- 
riched matter is assumed to mix homogeneously with the 
surrounding subcloud material (Fig. 2d). The remaining 
part is assumed to mix homogeneously either to a nearby 
subcloud hosting the next star formation event or to the 
ambient inactive cloud part (Fig. 2e). No mass-exchange 
is assumed between the subcloud and the ambient inac- 
tive cloud ISM during the time interval in which two or 
more star formation events occur within the same sub- 
cloud. In addition, subcloud material may be enriched by 
stars formed outside the subcloud. In this case, stars from 
elsewhere in the inactive cloud occasionally enter the sub- 
cloud region and enrich the subcloud by means of their 



ejecta (Fig. 2f). Stellar enrichment by old stellar genera- 
tions is assumed to proceed continuously with time but is 
considered in detail only at specific evolution times cor- 
responding to the occurence of any of the discontinuous 
processes referred to in Fig. 2. 

We define the subcloud core dispersal time Aidisp as 
the time between onset of star formation within a sub- 
cloud core region and the complete dispersal of this re- 
gion. This time interval constraints the mass of the most 
massive star that is able to enrich the subcloud core mate- 
rial before the core ultimately breaks up. Before dispersal 
of an entire subcloud, the subcloud will be enriched by 
the stellar populations it is hosting. After subcloud dis- 
persal (i.e. after a typical mixing time scale Ai m i x ), stars 
and gas belonging to the subcloud are assumed to mix in- 
stantaneously and homogeneously with the inactive cloud 
ISM. Subsequently, different cloud fragments may com- 
bine to form new subclouds wherein star formation occurs 
as soon as the critical conditions for star formation are 
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met. The mixing history of each subcloud determines the 
inhomogeneous chemical evolution of the inactive cloud 
part as well as that of nearby subclouds. For simplicity, 
we do not consider partial mixing of subcloud material to 
the inactive cloud. 

3. 2. Outline of model computations 

We perform Monte-Carlo simulations of the inhomoge- 
neous chemical evolution of a star forming gas cloud. The 
continuous process of formation and break up of subclouds 
and of the formation and dispersion of subcloud core re- 
gions associated with star formation, are followed as out- 
lined in the previous section. During the evolution of the 
cloud, we keep track of the total mass contained in gas and 
stars as well as the stellar and interstellar abundances of 
H, He, C, O, Fe, Mg, Al, and Si, both within each sub- 
cloud and the inactive cloud part. No instaneous recycling 
is assumed, i.e. metallicity dependent stellar lifetimes are 
taken into account. 

3. 3. Model input parameters 

Model input parameters for the reference model are listed 
in Table 1. We distinguish parameters related to: 1) the 
entire cloud and inactive cloud part, 2) active subcloud 
regions, and 3) individual star formation events: 

• Cloud and inactive cloud part: The initial cloud mass 
M c \ is treated as a mass scaling parameter (i.e. resulting 
abundances are not altered for different values of M c \). 
We here adopted M c \ = 5 10 10 M Q similar to that of the 
Galactic disk (e.g. Binney & Tremaine 1987). We assume 
a cloud evolution time t cv = 14 Gyr. This is comparable to 
the age of the Galaxy as derived from the age of the old- 
est globular clusters, i.e. 14±3 Gyr (e.g. Buonanno et al. 
1989). In our model, the impact of processes causing stel- 
lar abundance variations does not depend on the specific 
age of the Galactic disk assumed. 

We consider a total number of star formation events 
during the cloud evolution time of typically N s f = 100. In 
practice, N s { is limited only by the preferred model run 
time, i.e. 1-2 hours on a HP Apollo 715 machine. The total 
number of subclouds N sc \ is determined by the number 
of star formation events within each subcloud. For the 
reference model, we assume a maximum number of star 
formation events within one subcloud A s ™ ax = 1 so that 
N sc \ = N s f. Cloud initial abundances A qc i are as given in 
Table 1. Initially, the cloud is considered homogeneous, 
metal-free, and void of stars. 

• Active subclouds: In case of the reference model, we force 
subclouds to form at regular intervals of t cv / N sc \ = 1.4 
10 8 yr. We assume the subcloud mass M sc i directly pro- 
portional to the entire cloud gas-to-total mass-ratio \i at 
time of subcloud formation t sc \. This implies more massive 
subclouds to form at relatively high gas fractions fi(t). As- 
suming a constant star formation efficiency, this results in 



an exponential decrease of subcloud mass with disk age t 
(e.g. Clayton 1985): 

M scl = M scl (0) exp(-t/t dcc ) (1) 

where Ms C i(0) is the mass of a subcloud formed at t=0 
(which may vary between different models) and idec a 
characteristic time scale at which the mass of subsequent 
subclouds formed is assumed to decay (identical for all 
models). The assumption of a star formation rate (SFR) 
directly proportional to the subcloud formation rate is not 
essential for the results discussed here. 

The decay time scale idec is constrained observation- 
ally by the ratio of the average past to present SFR in 
the Galactic disk (<~ 3—7; e.g. Mayor & Martinet 1977; 
Dopita 1990). We here assume an exponentially decaying 
SFR with tdoc = 6 Gyr. As will be shown in Sect. 4.1, 
this SFR can account simultaneously for the actual gas- 
to-total mass-ratio in the disk of \i\ ~ 0.05—0.2 (Kulkarni 
& Heiles 1987; Binney & Tremaine 1987; see also Basu & 
Rana 1992), the smooth increase in the global AMR for el- 
ements such as O and Fe, and the magnitude of the current 
SFR in the Galactic disk (i.e. ~3.5 M Q yr -1 ; e.g. Dopita 
1987). In contrast, constant SFR models are inconsistent 
with these observations (Twarog 1980a; see also Clayton 
1985). 

The time between the formation and complete mix- 
ing of a subcloud to the inactive cloud part is defined as 
Ai m i x . This time scale has been considered to allow for 
the individual chemical evolution of a subcloud isolated 
from the inactive ISM (see below). 

• Individual star formation events: We assume the onset 
of star formation within each subcloud to coincide with 
the formation of the subcloud itself in case of the refer- 
ence model. This results in a grid of regularly spaced star 
formation times t s f ~ t sc \. 

We define the core dispersal time Aidisp as the time 
between onset of star formation t s f within a subcloud core 
and the moment star formation ends due to the actual 
break up of this core. Observational estimates of this time 
scale are generally ^10 7 yr (e.g. Garmany et al. 1982; 
Leisawitz 1985; Genzel & Stutzki 1989; Rizzo & Bajaja 
1994; Haikala 1995). For the reference model, we assume 
the entire subcloud to break up at time of dispersal of the 
star forming subcloud core, i.e. Ai m ; x = Aidisp- 

The star formation efficiency e> is defined as the 
amount of subcloud matter AM sc i turned into stars dur- 
ing star formation event j. In fact, the star formation ef- 
ficiency determines the amount of material to which the 
stellar ejecta of a previous stellar generation arc mixed 
within a given star forming cloud. Observational estimates 
for e in molecular clouds in the Galactic disk span a wide 
range: between a few tenths of a percent to ^50% (e.g. 
Wilking & Lada 1983). We will discuss the values assumed 
for e J in Sect. 4. 

Sequential stellar enrichment is taken into account by 
assuming that a fraction A-? of enriched material associated 
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with star formation event j is mixed homogeneously to the 
subcloud hosting the next star formation event. Subcloud 
material not converted into stars is mixed to the inactive 
cloud part after complete dispersal of the subcloud. The 
relative importance of these model parameters on the re- 
sulting stellar abundance variations will be discussed in 
Sect 4.2. 

3.4- Stellar evolution data 

We follow the stellar enrichment of the star forming cloud 
in terms of the characteristic element contributions of 
Asymptotic Giant Branch (AGB) stars, SNII, SNIa, and 
SNIb/c. This treatment is based on the specific abundance 
patterns observed within the ejecta of each of these stellar 
groups (e.g. Trimble 1991; Groenewegen & de Jong 1992; 
van den Hoek et al. 1996). We take into account metallic- 
ity dependent stellar element yields, remnant masses, and 
ages, while assuming the stellar ejecta to be returned at 
the end of the lifetime of the star (see e.g. Maeder 1992; 
Schaller et al. 1992). The respective time delays in enrich- 
ment by SNIa and SNII are accounted for in detail. A more 
detailed description of the combined set is given elsewhere 
(e.g. van den Hoek et al. 1996; see also electronic version 
of this paper). 

For AGB stars (initial mass m 8 M©) we adopt the 
metallicity dependent yields presented by Groenewegen & 
de Jong (1992). These yields are based on a synthetic evo- 
lution model for AGB stars and are succesful in explaining 
the observed abundances in carbon stars and planetary 
nebulae in the Galactic disk (Groenewegen et al. 1995). 
For Typc-II SNc, we use the explosive nucleo-synthesis 
yields (independent of initial metallicity) described in de- 
tail by Hashimoto et al. (1993) and Thielemann et al. 
(1993) for stars with 8^ m[M Q ] ^60. The 20 M Q model 
of this set accounts well for the observed abundances in 
SN1987A (Nomoto et al. 1991). Explosive nucleo-synthesis 
yields for Type-la SNe are adopted from Nomoto et al. 
(1984; model W7 for SNIa at Z = Z and Z = 0.0 of the 
accreted material; see also Yamaoka 1993) and for SNIb/c 
from Woosley et al. (1995). According to these yields, typ- 
ical amounts of iron produced are ~ 0.08 M Q for SNII, 
-0.8 M for SNIa, and -0.1 M for SNIb/c. 

The adopted yields for SNIa, b/c are relatively uncer- 
tain due to unknown details of the progenitor history and 
the explosion mechanism (either binary or single star evo- 
lution; see e.g. Smecker-Hane & Wyse 1992; Woosley et 
al. 1993). However, we do not believe that these uncer- 
tainties are relevant for the qualitative results obtained in 
this paper (cf. Sect. 5.2). 

Metallicity dependent stellar yields for stars during 
their wind (i.e. pre-SN) phase have been adopted from 
Maeder (1992, 1993), to whom we refer the reader also 
for a definition of the stellar element yields as used in the 
Appendix. For stars with m ^ 20 M we used the higher 
mass loss rates in case Z=0.02 (cf. Maeder 1992; Schaller 



et al. 1992). The mass m a of the helium core left at the 
end of the He-burning phase (or C-burning phase for mas- 
sive stars) has been used as input for the SNII and SNIb/c 
nucleosynthesis models referred to above. Yields were lin- 
early interpolated both in m and m a . Errors due to the 
coupling of these sets of stellar evolution data are proba- 
bly small and are neglected here (cf. van den Hoek et al. 
1996). Remnant masses and stellar lifetimes were adopted 
from the Geneva group as well (e.g. Schaller et al. 1992). 

For the reference model, the adopted IMF-slope, stellar 
mass limits at birth, and progenitor mass ranges for stars 
ending their lives as SNIa and SNII(+SNIb/c) are listed 
in Table 2. Stars with m > 60 M have been excluded 
because their theoretical yields are rather uncertain (e.g. 
Maeder 1992). We expect that the IMF- weighted contri- 
bution by such stars to the enrichment of the ISM is rel- 
atively low. Stars more massive than m^ Nn are assumed 
not to explode as supernova but to end as black hole (cf. 
Maeder 1992; Nomoto et al. 1994; Prantzos 1994; Tsu- 
jimoto et al. 1995). Consequently, stars with m <; m^ Nn 
contribute to the ISM enrichment during their stellar wind 
phase only. When no upper mass limit m^ Nn = 25—30 
M is introduced, models using up-to-date SNII yields 
predict abundances that are too high compared to those 
observed in the ISM, in particular for helium and oxygen 
(e.g. Twarog & Wheeler 1982; Maeder 1992, 1993; Timmes 
et al. 1995). 

We consider a fraction </> SNIa = 0.005 of all WD pro- 
genitors with initial mass between — 2.5 and 8 M to end 
as SNIa. A more detailed discussion of the contribution by 
SNIa is postponed to Sect. 5.1 (see also Ishimaru & Ari- 
moto 1995). In addition, we assume about one third of 
all supernova progenitors with mf Nn m^ Nn to end 

as SNIb/c, i.e. </, SNIb /c = [SNIb/c / (SNIb/c + SNII)] = 
0.33. This value is based on the observed ratio of current 
formation rates of SNII and SNIb/c in the Galaxy (e.g. 
van den Bergh & Tammann 1991; Tutukov, Yungelson & 
Iben 1992; Cappellaro et al. 1993). 

4. Results 

We present results for the inhomogeneous chemical evo- 
lution model described in the previous section. First, we 
consider the reference model which does not incorporate 
stellar abundance variations at a given age. Thereafter, 
we discuss models that do incorporate stellar abundance 
variations due to: 1) sequential stellar enrichment, 2) in- 
fall of metal-deficient matter, and 3) combined infall of 
metal-deficient matter and sequential enrichment. 

4-1. Reference model 

We consider a homogeneous gas cloud with initial condi- 
tions as listed in Tables 2 and 3. Within this cloud, active 
subclouds are formed at regular time intervals of 1.4 10 8 
yr so that in total N c i = 100 subclouds form during cloud 
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Table 1. List of input parameters (values listed for reference model) 



Cloud and inactive cloud part 



Mci 


5 10 iu M Q 


total cloud total mass 




14 Gyr 


cloud evolution time 


N st 


100 


total number of star formation events 


N scl 


100 


total number of subclouds formed 


X qcl (0) 


H=0.76 


initial cloud abundances; He=0.24 


For each subcloud i 


tscl 


1.4 10 s yr 


time at which subcloud is formed, i.e. each t cv / N sc \ — 1.4 10 s yr 


M scl 


exp 


exponentially decaying subcloud mass at time of formation (M sc \(t = 0) = 6 10 9 M Q ) 


jymax 


1 


maximum number of SF events within one subcloud 






mixing time of entire subcloud 


For each star formation event j 


t s C 


^scl 


evolution time at which SF-event j occurs 


Atdisp 


10 7 yr 


subcloud core dispersal time 


e 


0.50 


subcloud star formation efficiency 1 - 1 ' 


A 





efficiency of sequential enrichment^ 2 ' 



Notes: (1) e 3 is the mass fraction of the subcloud converted into stars during dipsersal 
time Ai^ isp , (2) \ j refers to the amount of stellar material returned to the subcloud 
core hosting the next star formation event. 



Table 2. IMF related parameters and stellar enrichment 



7 

mi, m u 

™SNII „,SNII 

mi , m u 

SNIa „SNIa 



m 



SNIa 
^SNIb/c 



m„ 



-2.35 slope of power-law IMF m 1 

0.1, 60 M stellar mass limits at birth 

8, 30 M© progenitor mass range for SNII and SNIb/c 

2.5, 8 Mq progenitor mass range for SNIa 

0.005 fraction of progenitors ending as SNIa 

0.33 fraction of SNII progenitors ending as SNIb/c 



evolution time t cv = 14 Gyr. We assume no time-delay 
between the formation of the subcloud and the actual on- 
set of star formation within that subcloud, i.e. t s { — t sc i, 
and further assume each subcloud to experience a single 
star formation event. During this event, lasting Atdisp = 
10 7 yr, half of the subcloud mass is converted into stars, 
i.e. e = 0.50. After each event, both gas and stars con- 
tained within the subcloud are mixed homogeneously to 
the inactive cloud part. We note that the assumption of 
e = 0.50 has no physical meaning here other than defin- 
ing the gas consumption rate as a function of cloud age. 
Model related quantitites are given in Table 3 (see Sect. 
5.2). 

Figure 3a shows resulting stellar and gas-to-total mass- 
ratios vs. age for the reference model. According to the 
assumed variation of subcloud mass with cloud evolution 
time (cf. Sect. 3.2), the gas-to-total mass-ratio decreases 
exponentially from fi c \ = 1 to 0.1 (corresponding decrease 
in subcloud mass converted into stars is shown in Fig. 3b). 
Figure 3b illustrates that the amount of gas returned by 



massive stars during each star formation event is less than 
~5% of the total amount of gas converted into stars during 
the same event. This ratio is determined primarily by t s { 
and e (see below). The amount of gas returned during 
Aidisp by the entire stellar population within the inactive 
cloud part vs. cloud age is plotted for comparison. 

Corresponding stellar and interstellar [Fe/H] and 
[O/H] abundance ratios are shown in Figs. 4a and b, re- 
spectively. At a given age, stellar and ISM abundances are 
exactly the same so that abundance inhomogeneities do 
not occur. Note that the resulting AMRs do not depend 
on the adopted value for M c i as long as the normalisa- 
tion of the SFR remains such that the condition of a cur- 
rent gas-to-total mass-ratio fj,\ of 0.1 is met. The reference 
model predicts [Fe/H] and [O/H] abundance ratios that 
are consistent with the mean EDV data for stars younger 
than ~ 10 Gyr. For stars older than 10 Gyr, agreement 
with the observations may be improved e.g. by considering 
cloud ages in excess of t cv = 14 Gyr or by detailed model- 
ing of the halo-disk enrichment at early epochs of Galaxy 
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Fig. 3. Reference model: a Stellar-to-total (dashed curve) and gas-to-total (solid line) mass-ratios vs. age. The gas-to-total 
mass-ratio for the inactive cloud part coincides with that for the entire cloud, b Subcloud mass AM sc i = eM sc i converted into 
stars (dashed curve) and amount of gas returned to the subcloud by newly formed stars within At s ! (solid curve) vs. age. We 
assumed M Bcl (t = 0) = 6 10 9 M©. Total mass returned by previously formed stellar populations present in the inactive cloud is 
shown for comparison (dash-dotted curve). Fluctuations in this curve result from integration over different time intervals, i.e. 
At B { and i sc i (cf. Table 1). 




-15 -10 -5 -15 -10 -5 -1 

Time [Gyr] Time [Gyr] [O/H] 



Fig. 4. Reference model: Stellar [Fe/H] and [O/H] abundance ratios vs. age. Stellar and ISM abundances at a given age are 
exactly the same, a [Fe/H]: stellar abundances at birth (full circles) coinciding with ISM abundances (solid curve). Each full 
circle represents a stellar generation with total mass of approximately AM sc i . Asterisks with error bars indicate the maximum 
stellar abundance variations observed in the Edvardsson et al. (1993a) data for main-sequence F and G dwarfs (averaged over 
1.5 Gyr bins). Data for stars older than 15 Gyr have been omitted because of incompleteness (see Sect. 2). b [O/H]: curves 
and data similar to those for [Fe/H]. c Theoretical [Fe/H] vs. [0/H]-relations including SNIa (full circles) and excluding SNIa 
(dashed curve). Mean [Fe/H] vs. [O/H] relation for the Edvardsson et al. data is shown as a straight line (horizontal marks 
indicate the observational range) 



evolution. We here concentrate on the stellar abundances 
observed during the last 10 Gyr of Galactic disk evolution. 

Our adopted values of </> SNIa = 0.005 and m^Nii = 30 
M Q provide optimal consistency with the mean observed 
[Fe/H] vs. [O/H] relation (cf. Fig. 4c). Clearly, the slope 
of the resulting [Fe/H] vs. [O/H] relation in case of enrich- 
ment by SNII only (m„ NI1 = 60 M Q ) is inconsistent with 
the observations. Thus, the data provided by Edvardsson 
et al. imply that SNII and SNIa,b/c nucleo-synthesis sites 
mixed their products together well. In addition, dilution 
of the supernova ejecta by more metal-deficient material 
is needed to comply with the range in [Fc/H] and [O/H] 
observed for F and G dwarfs in the SNBH. This is simply 



because theoretically predicted (lifetime-integrated) mean 
[Fe/H] ratios within the ejecta of supernova progenitors 
are in general much larger than those observed for long- 
living stars in the SNBH (see Sect. 3.4). Thus, whatever 
process is reponsible for the observed stellar abundance 
variations, both mixing of ejecta from different SN-types 
and dilution with mctal-dcficicnt material are involved. 

The resulting well-defined tight AMRs for the refer- 
ence model are similar to those predicted by conventional 
single-zone chemical evolution models (e.g. Twarog 1980a; 
Tinslcy 1980). Such models account for the global chem- 
ical enrichment of the Galactic disk ISM during the last 
10 Gyr, at least for elements like O and Fe, but they ob- 
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Table 3. Summary of model input parameters and resulting quantities related to the SFR (*scaled to M c \ = 2 10 M©) 



Model 


Fig. 


M 8C i(0) 


N sf 


£ max 


y m 


Mi 


SFRi 


INFi 


7t Nil 








r TV T 1 

[ M© ] 






r tv if i 

[ M© ] 




[ M© yi 


— li r i\ t 

' ] [ M© y 


— In r —In 
r ] [yr ] 


r —In 

[yr ] 


Reference 


3+4 


1.3 (10) 


100 


0.50 


30 


0.09 


2.9 




1.3 (-2) 


7.1 (-4) 


Sequential (single) 


5-1, 5-2 


2.1 (10) 


100 


0.95 


25 


0.11 


3.7 




1.5 (-2) 


8.2 (-4) 


Sequential (multiple) 


5-3 


1.5 (10) 


200 


0.50 


25 


0.21 


2.5 




1.0 (-2) 


6.4 (-4) 


Infall 


6 


6.3 (9) 


100 


0.95 


40 


0.16 


5.4 


3.1 


2.5 (-2) 


1.6 (-3) 


Infall+Sequential 


7-1, 7-2 


6.4 (9) 


136 


0.90 


25 


0.19 


5.2 


1.8 


2.1 (-2) 


1.3 (-3) 


Infall+Sequential 


7-3 


6.4 (9) 


166 


0.90 


25 


0.32 


4.4 


3.1 


1.8 (-2) 


1.1 (-3) 


Observations* 












0.05-0.2 


3.6±1. 


i>0.5 


2-4 (-2) 


;> 3 (-3) 



'References: 

M c i: Bahcall & Soneira (1980); Fich & Tremaine (1991). 

fii: Kulkarni & Heiles (1987); Binney & Tremaine (1987); see also Basu & Rana (1992). 
SFR: Dopita (1987); Walterbos (1988; based on IR observations); Mezger 1988 
INF: Mirabel & Morras (1984); Lepine & Duvert (1994). 

SNII & SNIa: van den Bergh & Tammann (1991); Tutukov et al. (1992); Cappellaro et al. (1993); Strom 
(1993). 



viously provide no explanation for the observed variations 
in stellar abundances at a given age. 

4-2. Sequential stellar enrichment 

In case of sequential star formation, efficient self- 
enrichment of a star forming gas cloud by successive stellar 
generations may result in abundance enhancements rela- 
tive to the abundances in the ambient ISM. When the 
local mixing time scale is larger than the time between 
two successive star formation events in such a cloud, these 
abundance enhancements can be deposited and recorded 
by newly formed stars. 

In our model, the impact of sequential enrichment on 
abundance inhomogencitics in the ISM is determined by: 
a) the dispersal time of the star forming region, b) the 
total number of stellar generations formed within one and 
the same cloud, c) the efficiency of sequential enrichment, 
i.e. the mass-ratio of the enriched stellar material and the 
cloud to which this material is mixed, d) details of stellar 
enrichment: e.g. the relative number of SNII and SNIa, 
and e) the IMF and stellar mass limits at birth. We dis- 
tinguish the effect of single and multiple sequential stellar 
enrichment on the stellar abundance variations. We will 
refer to single sequential enrichment as the case in which 
a star formation event induces subsequent star formation 
in a nearby cloud (when mixing enriched material to this 
cloud). 

4.2.1. Single sequential stellar enrichment 

We present results for models incorporating single sequen- 
tial stellar enrichment (N™ ax = 1). For illustration pur- 
poses, we consider only the alternating half of the sub- 
clouds to experience sequential enrichment. We define ini- 



tial masses of subclouds that are sequentially enriched as 
M sc \ = i)M* where M* is the total mass of gas converted 
into stars during the previous enriching star formation 
event. Initial masses of subclouds not involved with se- 
quential enrichment are assumed to decrease exponentially 
(as for the reference model). 

To maximize the effect of sequential enrichment on the 
stellar abundance variations, we assume $ = 0.2, e — 0.95, 
mf n = 25 M , (/> SNIa = 0.005, and an enrichment effi- 
ciency A = 0.95 (i.e. nearly all enriched stellar material 
ejected is mixed to the material wherein the next star 
formation event is induced). We consider cloud disper- 
sal times Atdisp ~ 10 7 yr. Such disperal times are among 
the largest ones deduced from observations of nearby star 
forming molecular clouds (see Sect. 3.3). Using a theoreti- 
cal age vs. turnoff-mass relation (e.g. Schaller et al. 1992), 
Atdisp can be related to the least massive star m enT able to 
enrich material before cloud dispersal. For instance, val- 
ues of A£ disp ~ 5 10 6 , 10 7 , and 2 10 7 yr, correspond to 
ra cm ~ 40, 15, and 12 M Q , respectively. 

Figures 5-1 and 5-2 (top and center panels in Fig. 5, re- 
spectively) illustrate the effect of sequential enrichment on 
the stellar abundance variations for cloud dispersal times 
of Aidisp = 10 7 and 2 10 7 yr, respectively. The extent 
to which sequential enrichment contributes to the stellar 
abundance variations is determined by the chemical evolu- 
tion of the ambient ISM. In general, large cloud dispersal 
times give rise to efficient locking up of metals in long liv- 
ing stars and enhanced stellar abundance variations rel- 
ative to the abundances in the ISM. Stellar abundance 
variations due to single sequential stellar enrichment are 
found to be maximal for Atdisp ^2 10 7 yr. Larger val- 
ues of A^disp allow stars less massive than m ~ 12 M Q 
to dilute the metal-rich ejecta of more massive stars (e.g. 
Hashimoto et al. 1993). 
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Fig. 5. Model results for single and multiple sequential stellar enrichment. Top panels: Single sequential enrichment assuming 
a cloud dispersal time tdisp = 10 7 yr. Center panels: Single sequential enrichment assuming fdi sp = 2 10 7 yr. Bottom panels: 
Multiple sequential enrichment (see text). Model results are shown for variations of stellar and interstellar [Fe/H] and [O/H] 
abundance ratios: a [Fe/H] vs. age, b [O/H] vs. age, c [Fe/H] vs. [O/H]. Stellar abundances are indicated by filled circles. 
Mean interstellar abundances (averaged over both active and inactive clouds) are indicated by solid curves. Average interstellar 
abundances within the inactive cloud ISM only (indicated by short dashed curve) do approximately coincide with the overall 
mean abundances. Remaining symbols and curves have the same meaning as in Fig. 4 



In case Aidisp = 2 10 7 yr, resulting stellar abun- 
dance variations due to sequential enrichment arc suffi- 
ciently large to explain the observed variations in [O/H]. 
In contrast, corresponding variations in [Fe/H] are much 
smaller than observed. This is true even though the mod- 
els presented here do account for sequential enrichment 
by SNIb/c which usually show theoretical [O/Fe] ratios 
much lower than SNII (e.g. Woosley et al. 1995. Results 
disagree with the observational fact that stellar abun- 
dance variations in [Fc/H] are considerably larger than in 
[O/H] (see also Gilmore & Wyse 1991; Edvardsson et al. 
1993a). An enhanced contribution of SNIb/c (i.e. assum- 
ing 0sNib/c > 0.33) or participation of SNIa to the process 
of sequential enrichment seems to be excluded by the ob- 



servations (van den Bergh & Tammann 1991; Tutukov et 
al. 1992). 

We conclude that single sequential stellar enrichment 
models are inconsistent with the observations because 
they: 1) result in [O/H] variations that are larger than 
those in [Fe/H], 2) predict current ISM abundances far 
below those observed, and 3) are difficult to reconcile with 
the apparent age independency of the stellar abundance 
variations (see Sect. 2). This conclusion is independent of 
the assumed cloud dispersion time scale Aidi S p> sequen- 
tial enrichment efficiency A, value of background level 
of iron-group elements set by SNIa (i.e. </> SNIa ), star for- 
mation efficiency e, and adopted IMF. Also, omitting the 
enrichment during the stellar wind phase of supernova pro- 
genitors does not alter this conclusion. 
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4.2.2. Multiple sequential stellar enrichment 

The effect of sequential stellar enrichment on the abun- 
dance variations among successive stellar generations can 
be very large, especially for high sequential enrichment ef- 
ficiencies and/or small amounts of cloud material to which 
the stellar ejecta are mixed before star formation is initi- 
ated. These conditions are naturally fullfilled when iso- 
lated gas clouds experience multiple star formation events 
(i.e. A r ™ ax >l) before mixing with the surrounding ISM. 
Since the gas content of an isolated cloud is reduced by 
each star formation event, cloud abundances rapidly in- 
crease when enriched by successive generations of massive 
stars. 

We consider two possible scenarios of multiple sequen- 
tial enrichment. In the first scenario, earlier generations of 
stars actually separate from the remaining subcloud mate- 
rial after dispersal of the subcloud core and do not further 
participate in the enrichment of the subcloud. Such mod- 
els result in substantial stellar abundance variations only 
under conditions and assumptions similar to those for the 
single sequential enrichment case. In the second scenario, 
all stellar generations formed in the subcloud continue to 
contribute to the enrichment until the entire subcloud has 
been dispersed. In this case, large stellar abundance vari- 
ations arise due to efficient recycling of the stellar ejecta 
from successive generations formed within the same cloud. 

We apply the second scenario and consider all sub- 
clouds to experience multiple sequential enrichment. As 
maximum number of star formation events within one and 
the same subcloud we assume iV™ ax = 4 (e.g. suggested by 
observations of OB subgroups in Orion; see Blaauw 1991). 
We adopt a sequential enrichment efficiency A = 1 (by def- 
inition within the same subcloud), a star formation effi- 
ciency e between 0.3 and 0.6 for each star formation event, 
and Atdisp = 10 7 yr. Other values of these parameters may 
provide similar results. Remaining model parameters are 
taken as for the single sequential enrichment model. 

In Fig. 5-3, we plot resulting stellar abundance vari- 
ations in case of multiple sequential stellar enrichment. 
Although variations caused by the first sequential en- 
richment event are relatively small (similar to the sin- 
gle sequential enrichment case assuming Atdisp = 10 7 yr; 
see Fig. 5-1), abundance variations caused by subsequent 
events can be as large as ~ 0.2—0.3 dex (depending on the 
subcloud abundances). As mentioned before, such large 
abundance variations are mainly due to ongoing sequential 
enrichment of the remaining cloud material by stellar gen- 
erations formed earlier in the cloud. We find that models 
incorporating multiple sequential stellar enrichment en- 
counter the same problems as single sequential enrichment 
models. However, the former models appear observation- 
ally far more justified. This is true in particular for the 
sequential enrichment and star formation efficiencies, as 
well as the cloud core dispersal times, required to obtain 
a given stellar abundance variation. 



We conclude that the observed stellar abundance vari- 
ations are difficult to explain by sequential stellar en- 
richment only. This conclusion is not altered when al- 
lowing for variations in sequential enrichment between 
distinct star formation events (e.g. by considering vari- 
ations in the IMF, relative formation rates of SNII and 
SNIb/c, etc.). Possible exceptions may be selective mix- 
ing of SNII nucleo-synthesis products to the material 
wherein star formation is induced, and/or cloud condi- 
tions that determine both the composition of the ejecta 
returned by a stellar generation (e.g. by means of the IMF, 
m^ Nn , Atdisp, contribution by SNIb/c, binaries) and the 
sequential enrichment efficiency (and/or star formation ef- 
ficiency). However, the impact of such conditions on the 
stellar abundance variations, which would imply that the 
IMF weighted SN yields used here would be considerably 
modified, is beyond the scope of this paper. 

4--S. Episodic infall of metal- deficient matter 

Infall of metal-poor material can account for abundances 
of newly formed stars which lie substantially below the 
abundances in the global disk ISM. In principle, abun- 
dance variations due to metal-deficient gas infall are de- 
termined by the abundances within the infalling gas, the 
gas infall rate, and the amount of disk ISM to which the 
infalling material is mixed before star formation is initi- 
ated. Stellar abundance variations due to infall of metal- 
rich material associated with SN ejecta from massive stars 
in the Galactic disk can be considered as a special case of 
sequential stellar enrichment and is not discussed here. 

Element abundances within the infalling gas are con- 
strained by the lowest abundances observed for disk stars 
with [Fc/H]<; —1. The Edvardsson et al. data imply infall 
abundances of [M/H] inf ^ -0.8 to -1.2 for e.g. M=C, Mg, 
Al, and Si. These abundances are consistent with observa- 
tions of interstellar clouds in the halo (see Sect. 5.2) and 
suggest that infall induced star formation is associated 
with the lowest abundances observed among disk F and 
G dwarfs in the SNBH. We assume infall abundances sim- 
ilar to the abundances observed among the oldest metal- 
poor disk stars, i.e. [Fe/H] = -1, [O/H] = -0.65, and 
a hydrogen mass fraction of X^0.72 (e.g. Bessell et al. 
1991). For simplicity, we do not account for abundance 
inhomogencitics within the infalling gas and assume infall 
abundances to be constant in time. Furthermore, we con- 
sider infall to occur as soon as the infall abundances are 
reached in the global disk ISM (presumably corresponding 
with the onset of star formation in the disk). 

The detailed manner in which gas infall varies with 
time is not essential for the results presented here as long 
as infalling gas plays an important role in determining the 
stellar abundances when it induces star formation. Here, 
we deal with the concept of infall induced star formation, 
i.e. the infalling gas initiates star formation when falling 
onto the disk. This concept is based on observations of 
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Fig. 6. Infall of metal- deficient material: model results are shown for variations of stellar and interstellar [Fe/H] and [O/H] 
abundance ratios: a [Fe/H] vs. age, b [O/H] vs. age, c [Fe/H] vs. [O/H]. Symbols and curves have the same meaning as in Fig. 
5 



infalling high-velocity clouds that strongly interact with 
disk ISM and initiate star formation therein as soon the 
critical density for star formation is reached (see Sect. 
5.2). Since star formation occurs by definition within ac- 
tive subclouds according to our model, infall is associated 
with subclouds experiencing star formation shortly after 
their formation. For simplicity, we assume each subcloud 
to contain an amount of infalling gas accumulated at the 
time star formation is induced. This amount is taken as 
a random fraction of the initial subcloud mass (i.e. be- 
tween and 1). In this manner, we allow for local and 
episodic gas infall onto the Galactic disk ISM. On aver- 
age, the gas infall rate is assumed to decay exponentially 
on a time scale tdec = 6 Gyr while its overall amplitude is 
constrained by m = 0.05-0.2 (cf. Eq. (1); see Table 3). 

We assume an initial disk mass M c \ = 3 10 10 M Q 
before the onset of gas infall. This results in a disk initial- 
to-final mass-ratio ( =0.55 according to an exponen- 
tional decrease of subcloud mass with cloud evolution 
time (M scl (0) = 2 10 9 M Q ; cf. Tabic 3). Clearly, stellar 
abundance variations due to metal-deficient gas infall are 
small at low levels of enrichment of the global disk ISM 
(i.e. large initial disk mass). Furthermore, continuous and 
large scale metal-deficient gas infall not associated with 
star formation results in relatively small stellar abundance 
variations. Thus, the magnitude of the stellar abundance 
variations is affected by the ratio of initial disk mass and 
total amount of infalling matter, as well as the amount of 
infalling gas that is involved with induced star formation 
in the disk ISM. 

Figure 6 displays resulting AMRs for iron and oxy- 
gen in case of episodic infall of metal-deficient gas. We as- 
sumed </> SNIa = 0.015, SNIb / c = 0.33, and mf 11 = 40 M . 
The value of m s u ml is taken larger than for the reference 
model (cf. Table 2) to obtain somewhat better agreement 
with the observations. Resulting abundances in the disk 
ISM follow the upper end of the abundances observed in F 
and G dwarfs younger than ~10 Gyr. Although the effect 
of global infall of metal deficient gas is generally to dilute 
the enrichment of the ISM, the inflow model results in 



larger ISM abundances than the reference model. This is 
mainly due to: 1) the asumption of a low initial disk mass 
which allows for a rapid early enrichment of the disk, and 
2) the assumption of local infall of metal-deficient gas and 
subsequent star formation therein so that infalling mate- 
rial has relatively small effect on the dilution of the global 
disk ISM. 

By varying the ratio of infalling matter and disk ISM 
within each subcloud, stellar abundance variations of A 
[Fe/H] ~ 0.8 dex and A [O/H] ~ 0.65 dex naturally can be 
accounted for. In addition, the scatter in [Fc/H] remains 
larger than in [O /H] since the iron abundances within the 
infalling gas relative to solar (i.e. [Fe/H]; n f = —1) are much 
smaller than that of oxygen ([0/H]; n f = —0.65). We note 
that the current gas infall rate of ~3.1 M Q yr" 1 predicted 
by the model shown in Fig. 6 is larger than suggested by 
the observations (see Sect. 5.1). However, other choices of 
model parameters, e.g. (, predict much lower infall rates 
while providing similar abundance results. 

Our models incorporating metal-deficient gas infall are 
in good agreement with the observed magnitude of stellar 
abundance variations and the slope of the [Fe/H] vs. [O/H] 
relation observed. However, these models predict current 
interstellar [Fe/H] and [O/H] abundance ratios of ~0.2 
dex above solar. This is in marked contrast with [O/H] 
abundance ratios of <~ 0.15 dex below solar observed both 
in interstellar gas and recently formed stars in the SNBH 
(see Sect. 2). In addition, these models appear to disagree 
with the observations on two other grounds. First, no sig- 
nificant scatter in the [Fe/H] vs. [O/H] relation is pre- 
dicted, contrary to what is observed for intermediate age 
disk stars (cf. Fig. 1). Part of the observed scatter may be 
due to experimental errors but variations of at least ±0.1 
dex in the [Fe/H] vs. [O/H] relation are probably real and 
have to be explained by any satisfactory model. A way 
to account for such scatter would be to allow for con- 
siderable abundance variations among different parcels of 
infalling gas. Secondly, these models predict stellar abun- 
dance variations to increase with time. This is inconsistent 
with the apparent constancy of the abundance scatter ob- 
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served. Possible ways out may be uncertainties in the ages 
of stars older than ~ 10 Gyr or disk evolution times in ex- 
cess of t ev ~ 14 Gyr. 

We conclude that models dealing with metal-deficient 
gas infall can probably be excluded as the complete an- 
swer to the stellar abundance variations observed, even 
though such models are in good agreement with both the 
observed abundance variations and [Fe/H] vs. [O/H] rela- 
tion. This conclusion is primarily based on the fact that 
such models predict mean current ISM abundances ~0.4 
dex larger than those observed in the SNBH. 

4-4- Metal-poor gas infall combined with sequential enrich- 
ment 

Motivated by the results previously discussed, we study 
the combined effect of sequential stellar enrichment and 
episodic infall of metal-deficient gas on the inhomogeneous 
chemical evolution of the Galactic disk. Such investigation 
is important also because these processes are observed to 
operate simultaneously in the SNBH (see Sect. 5.2). 

Attractive features of combined infall of metal-poor 
gas and sequential stellar enrichment are that a self- 
consistent explanation can be obtained for: 1) the presence 
of high metallicity stars at early epochs of star formation 
in the Galactic disk (due to sequential enrichment), 2) the 
presence of metal-poor stars at recent epochs of Galactic 
evolution (as a result of metal-deficient gas infall), 3) the 
nearly constant magnitude of the stellar abundance vari- 
ations during the lifetime of the disk, and 4) abundances 
in the local disk ISM that are currently below solar (as 
observed for oxygen). 

We show in Fig. 7 results for combined sequential stel- 
lar enrichment and metal-deficient gas infall. Model as- 
sumptions concerning each of these processes are similar 
to those described in the previous sections (e.g. e max = 
0.90, A = 0.95, mf n = 25 M Q , SNIa = 0.005, and 
0SNib/c = o 33. cf Tables 2 and 3). The three models 
shown in Fig. 7 differ only in the amounts of disk ISM in- 
volved with sequential stellar enrichment and infalling gas 
accreted during the lifetime of the disk. For each of these 
models, resulting stellar abundance variations and [Fe/H] 
vs. [O/H] relation are consistent with the observations. 
Clearly, models with combined sequential enrichment and 
metal-poor gas infall do not encounter the specific prob- 
lems involved when each of these processes is considered 
separately. 

We study the relative impact of metal-poor gas infall 
and sequential stellar enrichment on the resulting stellar 
abundance variations as well as the global enrichment of 
the ISM. First, we investigate the effect of varying the 
fraction of subclouds (i.e. the amount of star forming disk 
ISM) experiencing multiple sequential stellar enrichment 
(JV s 7 ax =4, Ai disp = 10 7 yr, e max = 0.9; see Sect. 4.2). This 
fraction increases from ~10% to ~25% when going from 
top to center models shown in Fig. 7. The remaining part 



of the subclouds is assumed to experience one single star 
formation event. Furthermore, we assume half of the sub- 
clouds to form stars partly from infall of metal-deficient 
gas (i.e. Figs. 7-1 and 7-2), regardless of the number of star 
formation events in each subcloud. Note that subclouds 
involved with metal-poor gas infall form predominantly 
stars with abundances below those in the global disk ISM. 
It can be seen that mean interstellar abundances and stel- 
lar abundance variations are not significantly altered when 
the fraction of ISM associated with sequential stellar en- 
richment is increased from 10 to 25%. However, when this 
fraction is further increased, more and more metals will 
be locked up in long living stars due to sequential enrich- 
ment and marked deviations from the observed [Fc/H] vs. 
[O/H] relation will occur (see Sect. 4.2). 

Secondly, we investigate the effect when the fraction 
of subclouds forming stars from metal-deficient gas infall 
is increased from 50 to 100%. In this case, stellar genera- 
tions are all formed according to infall induced sequential 
star formation and the total mass of infalling gas is in- 
creased by a factor two. This results in a reduction of the 
interstellar [Fe/H] and [O/H] abundance ratios by ~0.1 
dex (see Fig. 7-3). Interestingly, this marginally affects 
the magnitude of the resulting stellar abundance varia- 
tions (assuming e max = 0.9, A = 0.95) but strongly alters 
the number of stars with abundances below those present 
in the global disk ISM. We note that direct comparison of 
the abundance results with previous models is not justified 
because the enhanced gas infall model results in a current 
gas-to-total mass-ratio /Ui ^0.3, i.e. considerably higher 
than the fi\ = 0.1—0.2 indicated by the observations (cf. 
Table 3). To arrive at [i\ = 0.2, a reduction in initial 
disk mass from 3 10 10 to 2 10 10 M© would be required. 
In turn, this would result in ISM abundances and stellar 
abundance variations similar to that for models with more 
modest infall rates. 

4-5. Additional abundance constraints 

Keeping these results in mind, we study how models with 
combined sequential stellar enrichment and metal-poor 
gas infall behave when confronted with additional observa- 
tional constraints provided by the stellar abundance vari- 
ations and current ISM abundances of C, Mg, Al, and Si. 

Carbon abundance data for 85 F and G dwarfs in the 
SNBH have been presented by Andersson & Edvardsson 
(1994). These data show that there is a weak correla- 
tion between [C/H] and [O/H] (see Fig. 8). The shape 
of this correlation differs from that between e.g. [Fe/H] 
and [O/H]. In addition, the variation in [C/H] (i.e. ^0.6 
dex) at a given value of [O/H] is about three times larger 
than that in [Fc/H]. 

If the observed stellar abundance variations are caused 
by infall of metal-deficient gas only, one would expect that 
stellar abundances for all elements heavier than helium 
would be mutually correlated, e.g. similar to the correla- 
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Fig. 7. Model results for combined sequential stellar enrichment and episodic infall of metal-poor gas. Top panels: ~10% of 
the clouds is assumed to experience multiple sequential enrichment while half of the subclouds is involved with metal-deficient 
gas infall. Nearly 5% of the clouds undergoes both infall of metal-deficient material and sequential stellar enrichment. Center 
panels: 25% of the subclouds experiences multiple sequential enrichment while half of the subclouds is involved with metal-poor 
gas infall. Bottom panels: as center panels but all subclouds experience metal-deficient gas infall. Model results are shown for 
variations of stellar and interstellar [Fe/H] and [O/H] abundance ratios: a [Fe/H] vs. age, b [O/H] vs. age, c [Fe/H] vs. [O/H]. 
Symbols and curves have the same meaning as in Fig. 5 



tion between oxgyen and iron. In case of sequential stellar 
enrichment only, a similar behaviour would be expected 
only for elements that arc produced predominantly by 
SNII and SNIb/c. This implies that abundance-abundance 
variations between elements which are not synthesized 
predominantly within SNII and SNIb/c (such as C and 
N) , on the one hand, and elements that are produced pre- 
dominantly within supernovae (e.g. O, Si), on the other 
hand, may be conclusive about the importance of metal- 
deficient gas infall. 

We show in Fig. 8 results for: 1) infall of metal-deficient 
gas only (cf. Fig. 6), 2) multiple sequential stellar enrich- 
ment only (cf. Fig. 5-1), and 3) combined metal-poor gas 
infall and sequential enrichment (cf. Fig. 7-1). The in- 



fall model predicts no substantial scatter in the [C/H] vs. 
[O/H] relation but follows the trend in the observations 
well. Although the scatter in the [C/H] vs. [O/H] relation 
suggests that infall is not exclusively responsible for the 
observed variations, the shape of this relation indicates 
that infall is important. Conversely, the sequential enrich- 
ment model shows large scatter in the [C/H] vs. [O/H] 
relation but appears to deviate from the observed trend. 
The correlation predicted by the combined sequential en- 
richment and infall model appears in reasonable agree- 
ment with the observations. However, the observed carbon 
abundances at [C/H]~0. exhibit considerable more scatter 
than predicted by the model shown in Fig. 7-1 and seem to 
require a somewhat steeper increase of carbon relative to 
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oxygen. This may be due to variations in sequential stellar 
enrichment between different star formation events and/or 
variations in abundances within the infalling material. 

In our models, the overall shape of the stellar [C/H] 
vs. [O/H] relation is due to infall of metal-poor material 
with carbon abundances [C/0]i„f « —0.4. The reason why 
these relatively low carbon infall abundances are necessary 
to explain the observed trend in the [C/H] vs. [O/H] rela- 
tion, is unclear. A possible explanation may be a delayed 
carbon enrichment of the disk ISM, e.g. by low-mass SNIa 
progenitors that experience incomplete carbon burning or 
by low-mass AGB stars with small carbon yields. 

Sequential stellar enrichment seems inevitable to ex- 
plain the observed scatter in the [C/H] vs. [O/H] relation. 
Although large abundance inhomogeneities within the in- 
falling gas may reproduce the observed variations as well, 
such inhomogeneities in [C/O] appear inconsistent with 
the small scatter observed in e.g. [Fe/O]. Also, uncertain- 
ties in the derived carbon abundances may be consider- 
ably larger than those in O and Fe but are not likely to 
exceed ~0.2 dex (see Andersson & Edvardsson 1994; see 
also EDV). Therefore, it seems improbable that the scatter 
observed in the [C/H] vs [O/H] relation is due to observa- 
tional errors. Finally, it is difficult to see how chemical dif- 
ferentiation processes (e.g. dust depletion, element mixing 
to the surface, O/N-cycle, metallicity dependent nucleo- 
synthesis, etc.) can cause such large abundance-abundance 
variations among stars similar in mass and age. 

From the arguments above, we conclude that models 
incorporating both infall of metal-poor gas and sequential 
stellar enrichment provide an adequate explanation for the 
observed [C/H] vs. [O/H] relation. For the model shown 
in Fig. 7-1, we compare the predicted stellar abundance 
ratios Si, Mg, and Al vs. Fe with the observations in Fig. 
9. We find that slight offsets in [M/H] are present between 
the model predicted and observed relations. These offsets, 
most pronounced for Al, are due to details in the adopted 
stellar yields and related parameters (see Table 3), and are 
not essential for the following discussion. Interestingly, a 
number of observed features are naturally reproduced by 
these models. 

First, the observations suggest that variations in the 
stellar [M/H] abundance ratios decrease with increasing 
metallicity. This is theoretically predicted by the individ- 
ual effects of both sequential stellar enrichment and infall 
of metal-deficient material (as discussed above). For the 
model shown in Fig. 7-1, the resulting scatter in [M/H] 
at a given value of e.g. [Fe/H] is mainly due to sequen- 
tial enrichment (except at abundances [M/H]<^ —0.7) and 
is strongly related to the iron contribution by SNIa in re- 
gions that do not experience sequential enrichment. It can 
be seen that the scatter in [M/H] strongly decreases at so- 
lar metallicities and above. This may indicate that much of 
the observed element-to-element variations at high metal- 
licities is due to observational errors of ^0.1 dex in [M/H]. 
Alternatively, the observed scatter at high metallicities 



may imply that sequential enrichment by massive stars 
varies from one star formation event to another, e.g. by 
means of variations in the upper mass limit for SNII. 

Secondly, the predicted variation in [Mg/H] at a given 
[Fe/H] is substantially larger than that in [Si/H], consis- 
tent with the observations. In our model, this is due to: 
1) the fact that part of the Si comes from SNIa which are 
important contributors also to Fe (Mg is produced less 
efficiently in SNIa than is Si by about one order of magni- 
tude; e.g. Nomoto et al. 1984), and 2) the predicted ISM 
abundance of [Mg/H] is less by about 0.1 dex than that 
of [Si/H] (sec Fig. 9). 

Thirdly, the observed variation in [Al/H] is similar (or 
even larger) than that in [Mg/H] (see EDV). This be- 
haviour is also found for the model shown in Fig. 7-1. The 
predicted ISM abundance of Al is probably too low by 
~0.1 dex so that the resulting abundance scatter would 
be slightly reduced when correcting for this. The impact 
of sequential stellar enrichment is more pronounced for Al 
than for Mg and Si, due to the somewhat lower ISM abun- 
dance of Al (even after correction). Overall, we conclude 
that the magnitudes of the resulting stellar variations in 
Si, Mg, and Al vs. Fe appear in reasonable agreement with 
the observations. 

Comparison of variations in Mg, Si, and Al vs. O with 
the observations reveals a somewhat different picture from 
that vs. Fe. The observed variation of [Mg/H] with [O/H] 
has been shown in Fig. 1; Al and Si display a similar be- 
haviour. No trend is observed for variations in the scatter 
in the A[M/H] vs. [O/H] relation, in contrast to that in 
the A[M/H] vs. [Fe/H] relation. Furthermore, the observa- 
tions indicate mean variations in [M/Fe] of the same mag- 
nitude as those in [M/O]. In contrast, the model shown in 
Fig. 7-1 predicts variations in Mg, Al, and Si vs. Fe that 
are considerably larger than those vs. O. This is simply 
due to the fact that a substantial fraction of Fe originates 
from SNIa. Therefore, models predict hardly any scatter 
in e.g. [Mg/O] since both elements are synthesized pre- 
dominantly within SNII (and SNIb/c). This implies that 
either the scatter in the observed [Mg/H] vs. [O/H] rela- 
tion is due to observational errors or that an additional 
process is needed in the models to explain this scatter. 

In the former case, there would be no reason to believe 
the variations observed in the [M/H] vs. [Fe/H] relations 
either. However, we have argued above that at least part 
of this scatter is real. In the latter case, an additional 
mechanism causing the scatter in [M/O] could be vari- 
ations from one star formation site to another in the en- 
richment by SNII (and SNIb/c). Alternatively, such varia- 
tions could include variations in e.g. the IMF, upper mass 
limit for SNII, and/or the mass distribution of binaries. 
We have verified that such variations generally result in 
abundance-abundance scatter sufficiently large to account 
for the observed variations of 0.1 dex in [M/O] and suffi- 
ciently small to have a negligble effect on the scatter ob- 
served in [M/Fe]. Clearly, element-to-element variations in 



L.B. van den Hoek & T. de Jong: Inhomogeneous chemical evolution of the Galactic disk 



19 




Fig. 8. Comparison between observed and model-predicted [C/H] vs. [O/H] relation. Observations: data for F and G dwarfs 
in the SNBH from Andersson & Edvardsson (1995). Open circles represent stars with mean stellar galactocentric distances at 
birth within 0.5 kpc from the Sun (_Rq = 8.4 kpc). Full dots indicate stars with average distances within ~2 kpc from the Sun. 
Typical errors are indicated at the bottom right (top left panel). Model results: Predicted abundances of stars (full dots) and 
gas (solid line) for models incorporating metal-poor gas infall and/or sequential stellar enrichment 



enrichment from one star formation site to another would 
be a natural refinement of the sequential stellar enrich- 
ment models discussed before. 

We conclude that models incorporating both sequen- 
tial stellar enrichment and episodic infall of metal-poor 
gas provide a natural explanation for the observed stel- 
lar abundance variations and are consistent with the ISM 
abundances of C, O, Fe, Mg, Si, and Al. We find that 
the mean ISM abundances and abundance-abundance re- 
lations can provide only limited constraints on the relative 
importance of sequential enrichment and infall induced 
star formation in the Galactic disk. Therefore, improve- 
ments in observational and theoretical constraints are re- 
quired to disentangle the effects of these processes on the 
inhomogeneous chemical evolution of the Galactic disk in 
a more quantitative way. 

5. Discussion 

We briefly examine how the combined sequential stellar 
enrichment and metal-deficient gas infall models discussed 
in the previous section behave when confronted with inde- 



pendent constraints provided by the current star forma- 
tion rate in the Galactic disk and the chemical evolution 
of the Galactic halo. Thereafter, we discuss observational 
evidence in support of sequential enrichment and gas infall 
in the local disk ISM and consider possible implications of 
these processes for the chemical evolution of the Galaxy 
as a whole. 

5.1. Additional constraints 
5.1.1. SFR related constraints 

The combined sequential enrichment and metal-deficient 
gas infall (Fig. 7-1) predicts a present SFR of ^5.2 M Q 
yr _1 , and current rates of SNII (excluding SNIb/c) and 
SNIa of R sml = 2.1 1(T 2 yr" 1 and i? SNIa = 1.3 1CT 3 yr" 1 , 
respectively. These values are roughly consistent with the 
observations (i.e. within a factor of two; see Table 3). 
Adopted values of m^ NII =25 M and V SNIa =0.005 may 
be somewhat too low since the SN-rates scale with the 
predicted SFR. For the same model, the current gas in- 
fall rate is determined by the assumed disk initial-to-final 
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Fig. 9. Comparison between observed (left) and model-predicted (right) [M/H] vs. [Fe/H] relations: a) [Si/Fe], b) [Mg/Fe], 
and c) [Al/Fe]. Observations: data for F and G dwarfs in the SNBH from Edvardsson et al. (1993a). Symbols have the same 
meaning as in Fig. 8. Typical errors are indicated at the lower right of each panel. Model results: predicted abundances of stars 
(full dots) and gas (solid line) for model 7-1 



mass-ratio £=0.5 and by the time scale T; n f= 6 Gyr on 
which infall decays exponentially. This results in a cur- 
rent gas infall rate of 1.8 M Q yr _1 . Observations indicate 
a current gas infall rate of ~0.5 M yr -1 (e.g. Mirabel & 
Morras 1984). However, selection effects may account for 
an underestimate of a factor of 2—3 (sec Sect. 5.2). We 
note that higher values of £ and/or lower values of 7} n f 
may apply equally well. 



The predicted rates above all scale with the amplitude 
of the SFR. In turn, this amplitude is determined by the 
total cloud mass M c i =2 10 11 M Q and SFR decay time 
scale ideer ~ 6 Gyr assumed (see Sect. 3.2). Distinct val- 
ues of M c i and/or idecr will not affect the predicted stellar 
and interstellar abundances substantially, provided that 
a current gas-to-total mass-ratio n\ — 0.1 is maintained. 
We conclude that the adopted parameters for the model 
shown in Fig. 7-1 are consistent with observational con- 
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straints on the current SFR, gas infall rate, and supernova 
rates. Obviously, these observations do not yet provide 
tight constraints on e.g. idccr, Mi; an d C thus preventing a 
clear distinction between chemical evolution models based 
on these quantities (cf. Table 3). 

5.1.2. Constraints related to the enrichment of the Galac- 
tic halo 

The mean plateau value of [O/Fe]— 0.5 ±0.15 observed 
for halo stars with [Fe/H]^ —1 (e.g. Bessell et al. 1991) is 
presumably determined by the average [O/Fe] ratio within 
the ejecta of SNII (and SNIb/c) as well as the initial abun- 
dances within the halo ISM. For our models, the plateau 
value implies a maximum upper mass limit of SNII pro- 
genitors of m^ Nn «40 M Q , assuming initial mctallicities 
[Fe/H]< — 1, a Salpeter IMF, and stellar yields as de- 
scribed in Sect. 3.4. An even larger value for m^ Nn is im- 
plied when SNIb/c contributed substantially to the halo 
enrichment. 

We assumed m^ Nn = 25 M for the combined in- 
fall + sequential enrichment model (Fig. 7-1) discussed 
before. This results in [O/Fe]~0.2 at [Fe/H]^ -1 while 
omitting the contribution from SNIb/c would have re- 
sulted in [O/Fe] -0.25 at [Fc/H]£ -1. This is inconsistent 
with the observations. Possible solutions to this discrep- 
ancy are: 1) raf" and/or the IMF have changed between 
the time stars formed in the halo and the time of onset 
of star formation in the disk, 2) m^ Nn is actually ~40 
M Q for disk stars so that the predicted current disk ISM 
abundances of e.g. O and Fe increase and the effect of 
sequential enrichment is reduced, and/or 3) the adopted 
yields for SNII (and or SNIb/c) are in error at metallicities 
below [Fe/H]~ -1. 

Although the first two possibilities cannot be excluded, 
we favour the latter option since values of m^ Nn —30 M Q 
are suggested by recent models accounting for metallic- 
ity dependent yields of SNII in full detail (e.g. Timmes et 
al. 1995). We emphasize that the detailed yields for SNII 
at metallicities [Fe/H]^ —1 are not important for the se- 
quential enrichment and infall model results for disk stars 
presented in this paper but we just want to note here that 
it is difficult to explain the mean [O/Fe] ratio in halo stars 
using the same models. 

The observed breakpoint in the [O/Fe] vs. [Fe/H] re- 
lation at [Fe/H] 1.0±0.2 (e.g. King 1994) is generally 

associated with the time SNIa start to contaminate the 
global disk ISM (e.g. Gilmore & Wyse 1991; Bravo et 
al. 1993; Ishimaru & Arimoto 1995). In our models, the 
breakpoint in the [O/Fe] vs. [Fe/H] relation is mainly de- 
termined by: 1) the assumed fraction ^ SNIa = 0.005 of 
main-sequence stars with initial masses between 2.5 and 8 
M (e.g. Nomoto et al. 1984), 2) the delay time r SNIa = 
2.5 Gyr after which SNIa start to contribute to the enrich- 
ment of the ISM (e.g. Smecker-Hane & Wyse 1992; Ishi- 
maru & Arimoto 1995), and 3) the frequency distribution 



of SNIa as a function of age for a given stellar generation 
(assumed to be constant from r SNIa to r SNIa + 0.5 Gyr, 
and zero otherwise). These assumptions, in particular for 
r SNIa , strongly affect the increase of [Fe/H] with disk age 
and thus determine the scatter in and the slope of the 
[O/H] vs. [Fe/H] relation predicted. In fact, SNIa provide 
a background signal of iron-group elements on top of which 
stellar abundance variations due to sequential enrichment 
by SNII+SNIb/c occur. 

We assumed r = 2.5 Gyr for the models presented 
in this paper, as recently suggested by Ishimaru & Ari- 
moto (1995). However, this assumption implies different 
breakpoints in the [O/Fe] vs. [Fe/H] relation for models 
with disinct star formation and infall histories. The model 
shown in Fig. 7-1 predicts [Fc/H]=— 1 after —1.3 Gyr (and 
results [Fc/H] — —0.75 after 2.5 Gyr). This may indicate 
that the assumed value of r SNIa = 2.5 Gyr is in error. Es- 
timates for r SNIa based on stellar evolution calculations 
suffer from large uncertainties in the detailed evolution 
scenario for SNIa progenitors (see e.g. Smecker-Hane & 
Wyse 1992; King 1994) while theoretical estimates for 
r SNia based on the observed breakpoint may suffer from 
large uncertainties in the assumed chemical evolution of 
the halo (e.g. Ishimaru & Arimoto 1995). 

Alternatively, model assumptions related to: 1) the en- 
richment rate of the disk ISM by SNII (e.g. IMF and SFR, 
m^ Nn and SNII yields), or 2) the initial abundances of the 
material to which the SNII ejecta are mixed (e.g. disk mass 
at the onset of star formation, amount of gas infall, and 
infall abundances) may be in error. For instance, adopted 
iron yields for SNII at metallicities [Fe/H]sS —1 may be 
too high by about a factor of two. This would be consis- 
tent with the discrepancy in the [O/Fe] ratio for Galactic 
halo stars discussed above. Furthermore, large amounts of 
metal-poor gas infall would improve the consistency with 
the assumption of r SNIa = 2.5 Gyr. Other possibilities, 
such as higher values of the disk total-to-final mass-ratio 
or larger SFR decay times idocr, seem to be excluded by 
the observations (e.g. van den Hoek et al. 1996). 

We conclude that combined sequential enrichment and 
metal-poor gas infall models are consistent with the ob- 
served plateau value and breakpoint in the [O/Fe] vs. 
[Fe/H] relation provided that e.g. the adopted SNII yields 
at low metallicities [Fe/H]< —1 are too high by about a 
factor of two. At the same time, we have illustrated how 
sensitive our model results are to specific assumptions re- 
lated to the enrichment by e.g. SNII and SNIa. These as- 
sumptions may affect quantitive conclusions concerning 
the relative importance of sequential stellar enrichment 
and metal-poor gas infall. However, our qualitative con- 
clusion regarding the simultaneous presence of these pro- 
cesses in the local Galactic disk is not altered. 
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5.2. Observational support and implications for Galactic 
chemical evolution 

We briefly discuss observational evidence in support of 
sequential stellar enrichment and metal-deficient gas infall 
in the local disk ISM, and consider possible implications 
for the chemical evolution of the Galaxy as a whole. 

Sequential star formation has been argued to occur 
in nearby molecular cloud complexes including the well 
known Orion, Taurus- Auriga-Pcrscus, Cepheus, Carina, 
and Chameleon cloud complexes (e.g. see the review by 
Blaauw 1991; Megeath et al. 1995; Tcsti et al. 1995; Gold- 
smith 1995). In the Orion molecular cloud complex, age 
differences between OB subgroups are typically ^2—7 Myr 
(e.g. Genzel & Stutzki 1989; Cunha & Lambert 1992). 
When star formation proceeds on such time scales, mas- 
sive stars belonging to different OB subgroups may enrich 
the ambient molecular cloud material before a next round 
of star formation is initiated. Interestingly, the oldest OB 
subgroups in the Orion OBI association appear to have 
oxygen abundances that are lower by about 40% (~0.2 
dex in [O/H]) compared to the younger subgroups (see 
e.g. Olive & Schramm 1982; Cunha & Lambert 1992). 

Apart from observational support for sequential en- 
richment in nearby star forming regions, there are strong 
indications that the molecular cloud out of which the Sun 
formed has been enriched sequentially as well. Studies re- 
lated to extinct radioactive nuclides sue h as 53 Mn both 
in the Sun and in meteorites suggest that the protosolar 
molecular cloud has been enriched by high mass stars from 
a preceding OB association, about 10—25 Myr prior to the 
actual formation of the Sun (e.g. Cameron 1993; Swindle 
1993). 

Observations suggest that high-velocity inflow of 
metal-deficient gas towards the Galactic disk is a common 
phenomenon as well (e.g. de Boer & Savage 1984; Wakkcr 
1990; Schwarz et al. 1995). Estimates of the current gas 
infall rate onto the Galactic disk range from 0.2—0.5 M Q 
yr -1 , based on high- velocity clouds (HVCs; v J> 250 km 
s -1 ; e.g. Mirabel & Morras 1984; Lepine & Duvert 1994), 
to ^0.7 M.q yr -1 derived from the soft X-ray background 
(Cox & Smith 1976), and ~1.5 M Q yr -1 based on obser- 
vations of atomic hydrogen (Oort 1970). Since gas infall 
rates derived from the inflow of HVCs are likely underes- 
timates (e.g. Mirabel & Morras 1984), we estimate that 
~30— 40 % of the stars currently forming in the Galac- 
tic disk may be associated with infall (assuming a current 
SFR of - 3.5 M Q yr" 1 ; e.g. Dopita 1987). 

Observational support for infall induced star formation 
in the SNBH has recently been presented for the Orion 
cloud complex (Lepine & Duvert 1994; Meyer et al. 1994), 
the Gould Belt (e.g. Comeron & Torra 1994), and the ( 
Sculptoris open cluster (Edvardsson et al. 1995). These 
observations suggest that the most prominent star form- 
ing regions in the SNBH have been partly formed from 
infalling clouds from the Galactic halo. Circumstantial ev- 



idence for HVC impacts on the Galactic disk is based on: 
1) the existence of subgroups of young stars in star form- 
ing molecular clouds at high Galactic latitudes (like the 
Orion molecular cloud complex), 2) the displacement of 
OB star clusters with respect to the centers of their par- 
ent molecular clouds, 3) the alignment of the OB clusters 
in directions that are substantially inclined to the Galactic 
plane, 4) the age sequence of the aligned OB associations 
with an age of ~T0 7 yr for the oldest subgroups, and 5) the 
large elongated or filamentary structures observed in e.g. 
the Orion, Taurus, Monoceros molecular clouds connect- 
ing the clouds to the Galactic plane (see also Tenorio-Tagle 
et al. 1986; Franco et al. 1988; Gomez de Castro 1992). 
Many of these phenomena can be naturally explained by 
the interaction of a high velocity cloud with disk ISM and 
are difficult to reproduce by the process of sequential star 
formation (Lepine & Duvert 1994). Consequently, infall 
induced star formation appears to be a process frequently 
operating in the Galactic disk. 

We have argued that metal-deficient gas infall and se- 
quential stellar enrichment arc likely to play an impor- 
tant role for the inhomogeneous chemical evolution of the 
Galactic disk. These processes probably occur on scales 
ranging from single molecular clouds to the formation 
sites of open and globular clusters (e.g. Brown 1991). 
Large abundance variations of ~0.3 dex in [O/Fe] observed 
among metal-poor halo stars (e.g. Bessell et al. 1991; Nis- 
sen et al. 1994) suggest that sequential enrichment and gas 
accretion have been important for the chemical evolution 
of the Galactic halo as well. 

Nevertheless, combined metal-poor gas infall and se- 
quential stellar enrichment may be a too schematic pic- 
ture of the complex set of processes directing the chem- 
ical evolution of the Galactic disk. In particular, merger 
events with small companion galaxies may be important 
as well (e.g. Quinn et al. 1993). In such case, both gas 
and stars in the companion galaxy may add substantially 
to the observed abundance inhomogeneitics in the Galac- 
tic disk (Pilyugin & Edmunds 1995b). However, the sug- 
gestion that stellar populations from merging companion 
galaxies contribute substantially to the observed stellar 
abundance variations in the Galactic disk seems difficult 
to reconcile with: 1) the apparent homogeneous distribu- 
tion of these variations within the metallicity range ob- 
served at a given age of the disk, and 2) the small scatter 
observed in the element-to-element variations for stars in 
the SNBH. Notwithstanding, merging may be important 
for the chemical evolution of the Galaxy by adding large 
amounts of predominantly metal-poor material and initi- 
ating star formation in the disk. 

Inhomogeneous chemical evolution due to sequential 
stellar enrichment and/or metal-poor gas infall is probably 
important also in nearby galaxies such as the Magellanic 
Clouds and M31. In the Large Magellanic Cloud, large 
abundance variations of ^0.4—0.8 dex in [Fc/H] among 
similarly aged open clusters are observed (e.g. Cohen et 
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al. 1982; Da Costa 1991; Olsewski et al. 1991). Part of the 
variations may be accounted for by a radial gradient of 
-0.15 dex kpc- 1 in [Fe/H] (Kontizas et al. 1993). How- 
ever, the main part of these variations is likely due to 
triggered star formation in supershells as indicated by the 
close association of HII complexes with large HI holes ob- 
served in the LMC (Dopita 1985; Lortet & Testor 1988). 
In addition, gas infall may have affected the chemical evo- 
lution of the tidally interacting Magellanic Clouds. 

Observational evidence in support of shock-induced 
star formation by SNII in the spiral arms of M31 has been 
presented by Magnier et al. (1992). At these sites, young 
OB stars are observed to initiate recent star formation so 
that large abundance inhomogcncitics due to sequential 
stellar enrichment are expected, similar to those observed 
among OB associations in the Orion star forming cloud 
complex in our own Galaxy. 

Stellar and nebular abundance indicators reveal that 
substantial abundance fluctuations exist in the ISM of gas- 
rich galaxies (e.g. Roy & Kunth 1995). For instance, abun- 
dance inhomogcncitics in metal-poor galaxies such as IZw 
18 may be among the largest observed in external galaxies 
(e.g. Kunth et al. 1995) although this is still highly uncer- 
tain (Pettini & Lipman 1995). Whether the abundance 
fluctuations observed in dwarf galaxies are due to varia- 
tions in self-enrichment of the Hll-rcgions in these systems 
(e.g. Pilyugin 1992) and/or are related to selective loss 
of metals through galactic winds driven by massive stars 
(Roy & Kunth 1995; Martin 1996) is unclear. 

We expect that sequential stellar enrichment is gener- 
ally inefficient in dwarf galaxies because of their low gas 
densities, and that the effect of metal-poor gas infall on 
the stellar abundance variations is weak due to their low 
ISM abundances. Instead, star formation and abundance 
inhomogcncitics induced by metal-rich gas infall may be 
relatively important in these systems. 
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A. Appendix 

We describe the adopted model for the inhomogeneous chemi- 
cal evolution of a star forming gas cloud. The model can be 
applied to various mass scales, e.g. to the entire system of 
molecular cloud complexes in the Galactic disk or to the star 
forming core regions within a single molecular cloud. We start 
from a homogeneous, metal free gas cloud with a total mass 



M c i. At any evolution time t in its evolution, this cloud is sub- 
divided into N ac i active subclouds (with corresponding masses 
M* cl ) involved with star formation and an inactive cloud part 
(with mass M qc i) not involved with star formation. We assume 
matter to be freely exchanged within the inactive cloud part. 
Each subcloud i is formed at corresponding evolution times t l scl 
and is allowed to follow its individual star formation, mixing, 
and infall history. 

A.l. Model description, definitions and assumptions 

During the lifetime t cv of the star forming gas cloud a total 
number N s f star formation events is assumed occur. Each star 
formation event j presumably occurs within an active subcloud 
i. We define N* f as the total number of star formation events 
within subcloud i. For the reference model N° f = 1 and each 
star formation event j occurs in corresponding subcloud i = j. 
Subclouds are allowed to experience numerous star formation 
events, i.e. N s f > 1. During each star formation event j at time 
t = t J sf within subcloud i, a total mass of gas <5M* cl = e J M* cl (^ f ) 
is transformed into stars. 

We define At^ isp as the time between the onset of star for- 
mation within a subcloud core and the complete dispersal of 
this core region by supernova explosions and/or stellar winds. 
During At 3 disp the subcloud core is assumed to form stars. The 
profile of the star formation rate (SFR) during At diBp is as- 
sumed constant and identical for all star formation events. 
However, quantities such as the minimum stellar mass formed 
and IMF-slope are allowed to vary from one star formation 
event to another (cf. Sect. 4.2). The subcloud core dispersal 
time determines the mass of the most massive star that is able 
to enrich subcloud cloud material before the core ultimately 
breaks up. At time of core dipsersal, the newly formed genera- 
tion of stars has returned an amount of material <5M^ ot . Accord- 
ingly, the net amount of material converted into stars during 
star formation event j is given by: 5M^ f = e J <5M* cl (^ f ) — 5M^ et . 

Subclouds M* cl are formed from the inactive cloud ISM 
at cloud evolution times t = t z scl . When a subcloud forms it 
adopts the abundances of the inactive cloud ISM at t = t l scl . 
For each subcloud, we define a mixing time scale At^^ as the 
time between formation of the subcloud and the actual break 
up of the entire subcloud. The instant of break up of the sub- 
cloud may be either after one or more star formation events, 
or before star formation actually takes place. In this manner, 
material can be deposited within a subcloud region for a con- 
siderable period of time before being mixed to the surrounding 
ISM. The mixing history of each subcloud directs both the in- 
homogeneous chemical evolution of the inactive cloud and that 
of the neighboring subclouds. 

Before an entire subcloud breaks up its constituent mate- 
rial will be enriched by the stellar populations it is hosting. 
We assume the stellar enrichment of the subcloud to proceed 
homogeneously. In order to allow for sequential enrichment, 
we consider a fraction of enriched material ejected during 
star formation event j to mix homogeneously with subcloud 
core material hosting the next star formation event. Simulta- 
neous with the ejection of enriched material returned by newly 
formed stars, a substantial fraction of the ambient subcloud 
matter /t J M^ d may be swept up during dispersal of its star 
forming core. This subcloud material may mix to the subcloud 
hosting the next star formation event as well. The subcloud 
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hosting the next star formation event may be either the sub- 
cloud hosting the current star formation event or a subcloud 
nearby. No matter exchange is assumed between the subcloud 
and the surrounding ISM during the time between two star for- 
mation events occuring within one and the same subcloud. In 
case of the reference model, we do not consider mass transfer 
between subclouds, i.e. X 3 = k 3 = 0. 

After an entire subcloud breaks up its matter is assumed 
to mix homogeneously to the inactive cloud part. At the same 
time, stars associated with the dispersing subcloud become 
part of the stellar populations in the inactive cloud. After break 
up, different cloud fragments present in the ambient ISM may 
form new subclouds wherein star formation occurs as soon as 
the critical conditions for star formation are met. 

In addition to the individual chemical evolution of sub- 
clouds, which is directed by their star formation history and 
exchange history with the surrounding ISM, we allow for local 
enrichment of a given subcloud by stars that were not formed 
within that subcloud. This may be particularly important for 
low mass SNIa-progenitors which travelled considerable dis- 
tances from their birth sites and enrich their immediate sur- 
roundings at the time they explode as SNIa (cf. Fig. 2f; see 
below) . 

A. 2. Basic equations 

We keep track of the total mass of and abundances in stars 
and gas as a function of evolution time, both within each sub- 
cloud and the inactive cloud ISM. For each star formation event 
we use conventional chemical evolution model equations (e.g. 
Tinsley 1980, see below) except for including metallicity de- 
pendent stellar lifetimes, remnant masses and element yields 
(cf. van den Hoek et al. 1996). 

A. 2.1. Mass-exchange between subclouds and the inactive 
cloud ISM 

We denote AQ as the variation of a quantity Q between two 
cloud evolution times t — At and t. With M c i(£ = 0) the initial 
mass of the cloud and no stars initially present, i.e. M*(0) = 0, 
we can express the variations of mass of gas and stars within 
the cloud as: 

AM cl = Etf (t) AM' cl + AM qd (Al) 

AM, = zf'tW (AC{ - AEi) - AE«, qd (A2) 

where N sc i(t) is the current number of individual subclouds, 
N s f(i) the current number of star formation events within the 
cloud, AC the total mass of stars formed during star formation 
event j, and AE 3 the total mass of matter returned within 
At by stars formed during star formation event j. We recall 
conventional expressions for AC 3 and AE J (cf. Tinsley 1980): 

AC 3 = / mSj^Mjim) dmdt (A3) 

AE J = 

i-t rm (t-t e i-t diBp ) 

/ / (m - m rom (m)) Sj(t - r(m))Mj(m) dr^.M) 

J t — At Jm (t-t sf ) 

where Sj and Mj denote the SFR by number [yr _1 ] and IMF 
[M© -1 ] for star formation event j. For convenience, we ig- 
nored the index j for t B f, t diBp as well as for the stellar mass 



boundaries at birth m\, m u . We emphasize that both the stellar 
remnant masses m rom (m), lifetimes r(m), and turnoff- masses 
m (t) are a function of the initial metallicity Z, (containing 
all elements heavier than He) of the stellar generation under 
consideration. Variations in the total gas masses within the 
inactive cloud and subcloud i, i.e. M qc i and M* cl respectively, 
can be expressed as: 

AM q d = AE,, qd - E fo rmM s fc cl + EdispM^l (A5) 

AML, = [AE* - AC,] 1 + AM sf , prov - AM sf , ncxt (A6) 
AM scl = Ef = t (t) AMLi + E form M s fe cl - E disp MLi (A7) 

where AE„. qc i refers to the amount of material returned by 
stars present in the inactive cloud within time At. We followed 
both the stellar ejecta from recently formed stars within ac- 
tive subclouds and the ejecta from older stellar populations 
present in the inactive cloud ISM. The total amount of gas 
depleted by subclouds which are formed within time At is de- 
noted by Ef orm Mg cl . Similarly, the amount of gas returned by 
subclouds which become dispersed within time At is denoted 
by Edi sp Mg cl . We remark that the term between square brack- 
ets in Eq. (A6) refers to star formation events which occur 
within subcloud i. 

A. 2. 2. Supernovae Type la 

The term AE* in Eq. (A6) is related both to stellar generations 
which formed within subcloud i and to stars that entered the 
subcloud from elsewhere in the cloud. We will consider the case 
of SNIa progenitors stars only. Consequently, the term AE* can 
be expressed as two terms, i.e. AE* = E s f (AE j )* cl + AEsNia- 
The former term is related to star formation events which 
occured within subcloud i while the latter term is associ- 
ated with subcloud enrichment by SNIa-progenitors formed 
elsewhere in the cloud. We define the total amount of mat- 
ter returned by SNIa within subcloud i during time At as: 
AEsNia = ag NIa At_RsNiam rcm (m) where Rsma, is the total av- 
erage SNIa-rate in the entire cloud and a SNIa the corresponding 
fraction of SNIa that is is assumed to go off within subcloud i. 
In case of the reference model AEsNia =0. 

A. 2. 3. Mass-exchange between individual subclouds 

As matter may be transferred from one subcloud to another 
(or within one subcloud from one subcloud core to another) we 
include terms AM s ( iPrcv and AM s f jnoxt in Eq. (A6). The term 
AM B f iPre v corresponds to the amount of material added from 
the preceding star formation event to the core of the subcloud 
currently experiencing star formation. The term AM s f. next 
refers to the amount of matter mixed from the subcloud core 
actually experiencing star formation to the subcloud core host- 
ing the next star formation event. For each star formation event 
j which happens to occur in subcloud i within the time interval 
At we may write: 

AM sf , prcv = A^ 1 ^" 1 + K^Msd.prev (A8) 

AM sf , ncxt = A^M^+^M^i (A9) 

where M sc i iPrev is the mass of the subcloud hosting the pre- 
ceding star formation event. In this paper, we presented only 
results for k j = 0. In general, n 3 >0 has a similar effect as 
when reducing the sequential enrichment efficiency A^ . 
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A. 2.4. Chemical evolution of subclouds and inactive cloud 
ISM 

Expressions for the average abundance changes of element X 
within the entire cloud, inactive cloud part, and subclouds can 
be written as: 

A(X c iM c i) = Ef = t (t) A(X: cl ]V4i) + A(X qcl M qcl ) (A10) 

A(X qc lM qc l) = AEx.qcl — Eform-fqclMsd + EdispX^]M^.) 

A(X sc iM sc i) = 

T,^f t) A(X' Bcl Mld) + E form X qcl M s fc cl - E disp J^ ql M^A12) 
A(X s l cl M^) = [AE X - X scl AC«Y AM x , pro - AMx£4A3) 

where the meaning of each term can be found from its counter 
part in Eqs. A5-A7. Similarly, expressions for AMx, prC v and 
AMx.ncxt can be written as: 

AMx.prcv = A J '- 1 5M^- r 1 et + K J '- 1 (X BC iM sc i) pre v (A14) 
AMx.ncxt = AW Xret + K j X l sc iMi cl (A15) 

where AE^- is the total mass of enriched material of element X 
returned within time At by a stellar generation formed during 
star formation event j: 

*E? X = 

rt rm (t-t Bf -t diBp ) 

/ / AMx(m)Sj(( - T(m))Mj(m) dmdi (A16) 

Jt-At Jm (l-l sl ) 

AMx(m) — mpx(m) + {in {m))Xl (A17) 

and AMx {m) is the total mass of element X ejected by a star of 
initial mass m born with metallicity Xl during star formation 
event j. The term AMx(m) includes both newly synthesized 
stellar material and matter initially present at the time stars 
were formed. Initial stellar abundances Xl are determined by 
the abundances of the subcloud i (hosting star formation event 
j) at time t 3 Bf , i.e. Xi — X^ cl (t = £^ f ). Literature sources for the 
adopted theoretical metallicity dependent stellar yields px(w-), 
stellar lifetimes r(m), and remnant masses m re m(m) are given 
in Sect. 3.4. 
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